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ERRATUM 


All of the computations in the main report and in the appendix were 
performed with the following error in the programming of the Newton- 
Raphson solution of the burning rate equation: PONOZ was used in place 

of PNOZ in Table IV-2, p. 43, line l6 of the report. The program statement 
is correct in the report, but because of the error numerical results 
obtained from the program will differ slightly from those presented in the 
report expecially during the initial portion of the traces. The variations 
are not large enough to change the conclusions of the report or to invalidate 
the use of the design charts presented in the appendix. 


ACKNOWLEDGEMENTS 


The author expresses grateful appreciation to 
Mr. J. M. Lyon, Graduate Research Assistant in 
Aerospace Engineering, Auburn University, who skill- 
fully and enthusiastically programmed the analysis 
given in this report for the digital computer and 
rendered other valuable assistance. 

The advice and helpful suggestions given by 
personnel at the George C. Marshall Space Flight 
Center who used the computer program during its 
development and provided data for analysis are also 
acknowledged. The aid thus received from Messrs. 

B. W. Shackelford, Jr. (NASA Project Coordinator) 
and K. W. Jones is especially appreciated. 


ii 



TABLE OF CONTENTS 


Acknowledgements . . . . ii 

List of Figures iv 

List of Tables vi 

Nomenclature vii 

I. Introduction and Summary 1 

II. Discussion of Input and Output k 

III. Analysis ....... 19 

IV. The Computer Program 36 

V. Test Cases j6 

VI. Concluding Remarks .... ...... 89 

References. ........ 90 

Appendix 91 


iii 



LIST OF FIGURES 


Figure II-l. Basic motor dimensions 6 

Figure II-2. Standard star grain cross-section . 11 

Figure II-3. Truncated (slotted tube) star grain cross-section ... 12 

Figure IV— 1. Flowchart for main program . 60 

Figure IV-2. Flowchart for area subroutine ..... 68 

Figure IV-3 . Flowchart for output subroutine ... 75 

Figure V-l. Sketch of SKM Number 1 . 79 

Figure V-2. Pressures and vacuum thrust versus time for 

SRM number 1 82 

Figure V-3. Sketch of SRM number 2 . 8H 

Figure W-h. Head end pressure and vacuum thrust versus time 

for SRM number 2 . . . . . . . . . 87 

Figure A-l. Nozzle end stagnation pressure versus time for 

various nozzle radial erosion rates (SRM l) 9U 

Figure A-2. Nozzle end stagnation pressure versus time for 

various initial grain temperatures (SRM l) 95 

Figure A-3. Nozzle end stagnation pressure versus time for 

various numbers of burning slot faces (SRM l) 96 

Figure A-1+. Nozzle end stagnation pressure versus time for 

various numbers of star points (SRM l) 97 

Figure A-5. Nozzle end stagnation pressure versus time for 

various lengths of star grain (SRM l) 98 

Figure A-6. Nozzle end stagnation pressure versus time for 

various differences in web thickness x^, over the aft 

end tapered grain length (SRM l). . . . a 99 

(Continued) 

iv 



Figure A- r 

Figure A-l 
Figure A-| 


LIST OF FIGURES (Cont'd) 


. Nozzle end stagnation pressure versus time for 

various differences in Veb thickness z 0 between the 
ends of the main grain (SRM l) 

. Nozzle end stagnation pressure versus time with 

and without thrust termination ports (SRM l) ..... . 

. Nozzle end stagnation pressure versus time for 

various values of the erosive burning constants (SRM l) 


. 100 

. 101 

. 102 


v 



LIST OF TABLES 


TAVLE IV-1. Data Card Formats 38 

TABLE IV-2 . SRM Design and Performance Analysis 1+0 

TABLE V-l . Input Data for SRM No. 1 ... 80 

TABLE V-2 . Sample Printout of Transient Values for SEM No. 1 . . . 8l 

TABLE V-3. Printout of Weights and Related Values . . . . . . . . .83 

TABLE V-4. Input Data for SRM No. 2 85 

TABLE V-5. Sample Printout of Transient Values for SRM No. 2 ... 86 

TABLE V-6. Comparison of Computed Weights and Related 

Values with Results of Others ... ...... 88 


vi 



NOMENCLATURE 


An asterisk before the symbols indicates an input variable. All input 
subscripted and non-subscripted variables are listed separately. 


English 


Symbol 

Definition 

Units Used 

*a 

Propellant burning rate coefficient 

in/sec-psi n 

A 

Cross-sectional or flow area 

in 2 

■\>m 

Main burning surface evaluated during tailoff 

in 2 

■\in ,A bp » 

Total burning surface associated with nozzle 

in 2 

#A bs 

end, port and slot surfaces, respectively 


*Vr 

Burning surfaces (as a function of y) at nozzle 
end for tabular input 

in 2 

* /l bpT 

Burning surfaces (as a function of y) of port 
sides for tabular input 

in 2 

*(y 

H3ST 

Burning surfaces (as a function of y) of slots 
for tabular input 

in 2 

*\to 

Sliver burning surface evaluated during tailoff 

in 2 

A Gs 

Initial burning area of star grain flat ends 

in 2 

A 

P 

Port area 

in 2 

*A 

phT* pnT 

Controlling port areas at head and aft ends, 
respectively (as a function of y) for purely 

in 2 


tabular inputs 

in 

*B ,*B 
bp bs 

Additive burning surfaces input as special 
equations for port, slot and nozzle surfaces, 
respectively 

in 2 

*c* 

Characteristic exhaust velocity 

ft/sec 

C F 

Thrust coefficient 
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NOMENCLATURE (Continued) 


English 

Symbol Definition Units Used 


D 


Diameter 


in. 


D* ,*D. * 
1 


Instantaneous and initial diameters, 
respectively, of the nozzle throat 


in . 


*D 

cc 


Estimated mean diameter of the cylindrical in. 

portion of the case 


*D 

e 


Diameter of the nozzle exit 


in. 



Length average initial diameter of controlling in. 

length of circular perforated grain 


*D 

o 


Length average outside diameter of circular in. 

perforated grain, excluding lengths in closures 


D^,p Initial diameter of the thrust termination in. 

passageways in the grain 


*f 

*E 


F 

g 


o 



h 


#h b 

*h 

c 

*1 


I 

J 


op 
sp * 


Fillet radius at star valley 

Radial erosion rate of the nozzle throat 

Thrust 

Gravitational constant 

Integer designating type of grain perforation 

Instantaneous vehicle altitude 

Estimated burnout altitude 

Axial length of nozzle closure 

Integer designating type of program input 

Specific and total impulse, respectively 

Throat-to-port cross-sectional area ratio 


in. 

in/sec 
lbf . 
ft/sec 2 


ft. 


ft. 


in. 


lbf-sec/lbm,lbf-sec 
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NOMENCLATURE (Continued) 


English 


Symbol 

Definition 

Units Used 

CM 

t^ H 

* 

Empirical constants in nozzle weight 
estimation equation (Eq. 63) 

lbm/in 3 , 

in/sec 

* K eh- K en 

Erosion rate of insulation at head and 
nozzle end of case, respectively 

in/sec 


Distance between center line of motor and 
fillet center of standard star grain 

in. 

l 

Vi 

Length of sides of truncated star point 
excluding fillet 

in. 

*0 

TP 

Initial length of termination passages between 
centers of gravity of perimeters of bases 

in. 

*L 

Total initial perforated grain length including 
gaps at slots 

in. 

*L 

cc 

Length of cylindrical portion of case 
including forward and aft segments 

in. 

L Gc’ L *Gci 

Instantaneous and initial total axial lengths, 
respectively, of circular perforated grain 
(not including gaps) 

in ♦ 

* L Gni 

Initial axial length 

in. 

* L Gsi 

Initial total length of star-shaped perforated 
grain 

in. 

#L Ta 

Estimated length of grain at the aft end at start 
of first tailoff having an additional taper 
not represented by z q or 0 G 

in. 

m 

Mass 

slugs 

M 

Mach number 


*n 

Burning rate exponent 


*n lS % 

Number of standard derivations in maximum 
chamber pressure and case yield strength, 
respectively, to be used as bases for design 
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NOMENCLATURE (Continued) 


English 

Symbol 


*n 


n 


*n 


*n 


n seg 

% 

P 

*Q 

r 

R 

*R 


op 


eL 


*R 


*S 


G 

*S 


op 


y 

*s 


yc nom 


*S 


yn 


Definition Units Used 

Number of burning flat end surfaces of a star 

grain located at the extreme nozzle end of 
the chamber 


Number of star points 

Number of burning flat end surfaces of a star 
grain not located at the nozzle end of the chamber 


Number of case segments 

Number of thrust termination ports 

Pressure lbf/in 2 

Integer designating grain arrangement 

Burning rate in/sec 

Radius in . 

Outside radius of a star grain in. 


Initial Reynolds number of gases in propellant 
port based on port length 


Initial radius of truncated star points 

Number of burning flat slot sides of a 
circular perforated grain not including the 
nozzle end 

Burning perimeter of a star grain 
Integer designating type of star grain 
Yield strength 

Nominal yield strengths of the case and 
nozzle structural materials, respectively 


in. 


in. 


lbf/in 

lbf/in 


x 



NOMENCLATURE (Continued) 


English 


Symbol 

Definition 

Units Used 

t 

Time 

sec. 

t b 

Calculated total operating time of motor 

sec . 

%1 

Estimated time at burnout 

sec . 

S 

Time at which propellant burning ceases 

sec. 

T £ 

Loss in length of burning propellant associated 
with additional tailoff controlling taper near 
the nozzle end of the grain 

in. 

u 

Velocity 

ft/sec 

V ,*V ., 

C Cl 

*v 

cf 

Instantaneous, initial and final volume, 
respectively, of chamber gases 

in. 

w 

Weight 

Xbm. 

W . ,W _ , 
pi* p2 * 

W 

P 

Weight of propellant burned based on mass 

discharge rate, volume calculations, and 

the arithmetic average of W and W _ 

pi p2 

lbm. 

X 

out 

Distance burned at which propellant breaks up 

in. 

Ta 

Difference in web thicknesses at ends of 1^ 

in. 

X 

Ratio of the sum of the grain burning areas 
excluding nozzle end burning areas to the 
total sum of the burning areas 


y 

Distance propellant has burned from initial 
surface 

in. 

*z ,z 
o’ 

Initial and instantaneous differences , 
respectively, between web thickness at head 
and nozzle ends of controlling grain length , 
excluding any initial length associated with 

Sa ° r 6 G 

in. 
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NOMENCLATURE (Continued) 


Greek 

Symbol 


'a 


eb 


n 


*£3 

*Y 

r 


* 6 c’*V 


ins 


&A bT 


# &D i 


&P c /&y 


* ( AP /Ay ) 


AT. 

i 


Ay, Ay, 


out 


Definition 

Erosive burning coefficient in the 
Robillard-Lenoir rule 

Nozzle exit half angle 


Erosive burning pressure coefficient in the 
Robillard-Lenoir rule 

Specific heat ratio of propellant gases 

A function of y given by Eq. k 

Specific weights of case, liner and case 
insulation, respectively 

Burning surface area change resulting from thrust 
termination passageways 

Difference between initial circular perforated 
grain diameter at nozzle end of controlling 
length and D^ 

Rate of change of chamber pressure with respect 
to distance burned 

Depressurization rate at which propellant is 
extinguished 

Maximum expected increase in initial grain 
temperature above design point condition 

Incremental distance burned and initial value of 
same 


Units Used 


in 2 *®-ft° ^/sec 1 , 8 lbf° 


radians 


lbm/in 3 




xn. 


lbf/in 3 

lbf/in 3 


in. 


e 

«F 

5 M 

*e ,*e , 

cn ch 


Nozzle expansion ratio 

Thrust loss coefficient 

Motor mass fraction: i,e. , W /W.. 

’ p M 

Approximate acute angle bonded circular 
perforated grain makes with motor center line 
at the head and nozzle closures, respectively. 
Also referred to as the grain contraction angles , 


radians 



KOMEWCLATURE (Continued) 


Symbol 

Definition 

Units Used 

* 9 f 

Angular location of fillet center of 
standard star from line of symmetry 

radians 

* 

a> 

O 

Angle burning surface element of circular 
grain makes with longitudinal axis of motor 
at the nozzle end of the chamber at the nozzle 

radians 

* 9 n 

Hozzle cant angle 

radians 

P< 

CD 

* 

Angle of standard star point 

radians 

e 

& 

Half angle of star grain sector 

radians 

6 ■. 
si 

Angular location of slot side of truncated 
star grain 

radians 

# fi 

8 tp 

Acute angle between axis of thrust 
termination passage and motor axis 

radians 

*» 

Absolute viscosity of chamber gases 

Ibf-sec/in- 

•'"p’k 

Temperature sensitivity coefficient of 
pressure at constant K = A^/A* 

/°F 

% 

Solid propellant density 

slugs /in 3 

P S yc yc 

Statistical variation in P and S y respectively; 
i.e. } one standard deviation/the nominal value 
of the variable 


% 

Thickness of the liner 

in* 

♦t 

s 

Thickness of propellant web at slot bottom of 
truncated star grain 

in* 

#T 

Teff 

Estimated "effective" web thickness of 
termination port 

in. 

*T 

W 

Web thickness of main propellant grain 

in. 

*T 

VS 

Web thickness of standard star (same as r w 
except for some combination grains) 

in* 


xiii 



HOMES CLATUKB (Continued) 


Greek 



a 


c 

d 

D 

e 

f 

G 

h 

i 

ins 

«j 

a 

M 

max 


Definition Units Used 

Safety factor on nozzle ablative material 

Safety factor on case thickness 

Safety factor on nozzle insulation thickness 

Safety factor on nozzle structural material 

Additional or ambient 

Arithmetic averages between head and nozzle 
end values 

Case, chamber value or cylinder 
Dome 

Discharged 

Exit or erosion 

Final 

Generated 

Head end of grain 

Initial 

Insulation 

Incremental index; i.e., number of increments burned 

Liner 

Motor 

Maximum value 

xiv 



HOMENCLATUBE ( Concluded 


Subscripts 

meop 

n 

nc 

ns 

nT 

nom 

o 

pc 

ps 

pa? 

sc 

ss 


( Cont 1 d ) 

Maximum expected operating condition 

Nozzle or nozzle end of grain 

Nozzle end of circular perforated grain 

Nozzle end of star grain 

Nozzle end tabular value 

Nominal value 

Stagnation 

Port of circular perforated grain 

Port of star grain 

Port tabular value 

Slot of circular perforated grain 

Slot of star grain 


sT 

vac 

v 

Superscripts 

* 


Slot tabular value 
Vacuum 

Value at veb time 
Choked throat value 


Arithmetic average value over Ay 
Time rate of change 


xv 



I. INTRODUCTION AND SUMMARY 


This report presents the results of research performed at Auburn 
University during the period from June 1 6 to October 12 , 1972, under Modi- 
fication No. 6 to the Cooperative Agreement, dated October 6, 1969, between 
the NASA Marshall Space Flight Center and Auburn University. The principal 
objective of the research was to develop a simplified computer program 
which would serve as an aid in preliminary design of solid-propellant 
rocket motors (SRMs). The program would also permit generation of SRM 

parametric and performance data which were to be compiled as a part of 
the task. 

A number of computer programs are already available to government and 
industry for use in SRM design which will give results based on quite precise 
mathematical formulation of burning surface geometry. The general tendency 
has been toward developing more and more sophistication in these programs. 

This creates substantial difficulty for the user in understanding the program 
and preparing inputs; lengthy orientation times are frequently required. Also 
use of such programs involves long operating and turn-around times on the 
computer and attendant high costs . Elementary textbook -type approaches are 
of course possible, but the accuracy requirements must be considered. Indeed, 
in view of the high cost associated with SRM development, a propulsion 
contractor would be ill-advised in not making some use of a tried and proven 
sophisticated program available to him. On the other hand, it is questionable 
whether the accuracy obtained with the sophisticated programs is sufficient 
to justify their adoption for the bulk of preliminary design work. Unknowns 
such as the erosive burning characteristics of a new propellant, or even a 
new SRM design with a well characterized propellant, cast doubt on the 
accuracy of performance predictions. In view of this, a less precise approach 
is often acceptable in preliminary design; for example, in determining the 
effects to anticipate from a given design change. 


This report presents an analysis and a computer program which represents 
a compromise between the more sophisticated programs using precise burning 
geometric relations and the textbook-type of solutions. The program requires 
approximately 900 computer cards including a set of 20 input data cards required 
for a typical problem. The computer operating time for a single configuration 
is approximately 1 minute and 30 seconds on the IBM 360 computer. About 1 
minute and 15 seconds of the time is compilation time so that additional 
configurations input at the same time require approximately 15 seconds each . 

The program uses approximately 11,000 words on the IBM 360. The program is 
written in FORTRAN IV and is readily adaptable for use on a number of different 
computers; e.g., IBM 704b, IBM 709 1 * and Univac 1108. 

Based on input of solid— propellant rocket initial geometric properties , 
data on the characteristics of the propellant and inert materials, and design 
criteria constants, the program computes the weights of major components and 
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performance characteristics of the SRM including specific impulse, thrust, 
chamber pressures, and other significant time varying parameters. The program 
may be used for rocket motors of all sizes; however, primary consideration 
was given in this study to very large SRMs (120 in. dia. and above). There- 
fore, some of the empirical expressions for component weights developed, such 
as the one for nozzle weights, may be more accurate for larger motors than 
for smaller ones as discussed further in the body of the report. 

With large motors again in mind, particularly those under consideration 
for application to the Space Shuttle Booster (References 1-li), the program 
was designed to treat segmented configurations with both star and circular 
perforated grains present in various arrangements in the same SRM. In 
addition, monolithic grains and either standard star or truncated star (slotted 
tube) shapes may be analyzed. While there is considerable flexibility 
provided in specifying the geometry of the ends of circular perforated grains, 
there are severe constraints on the specification of the geometry of star 
grains and on passageways in the grain for thrust termination ports. To 
offset this, as a program option, the user may specify tabular data to 
correct precisely or approximately for variation in burning surface geometry 
with distance burned. Alternately, the entire burning area geometry may be 
specified by tabular input. 

The key to the relative simplicity of the program, aside from the use 
of end constraints on burning surface geometry, is the way in which the 
variation in surface regression along the length of the grain is treated. 

The more sophisticated programs result from mathematically dividing the 
grain axially into a large number of segments. For each segment, in turn, 
the governing fluid dynamic equations are solved. The difficulty here is 
that the flow velocity at each segment is unknown until the burning rate is 
determined, but the burning rate depends on the velocity. Therefore, an 
iterative technique must be used which requires assumption of the conditions 
at the head end of the grain, solution of the governing equations for each 
segment in turn and finally a comparison of the total amount of mass generated 
with that theoretically discharged through the nozzle. If the check fails, 
a new assumption and repetition is required to obtain a solution for one 
instant of time. The advantage of this method is that it permits the surface 
regression to be determined at each segment for each instant of time, but 
it is clear that this results in both long computation time and large 
computer storage capacity . The present analysis avoids the above mentioned 
difficulties through the assumption that the surface regresses in a linear 
fashion along the length of grain consistent with the difference in surface 
regression at the head and aft ends for which regression (burning) rates are 
separately calculated. Support for this assumption as an approximation may 
be found in the results of interrupted burning tests (See for example Reference 
5, p. 1+18) . As already noted, the nature of erosive burning characteristics 
is at best poorly known which casts some doubt on a need for a more rigorous 
model of the burning surface geometry. Separate determination of head and 
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nozzle end burning rates taking into account erosive burning effects provides 
a means of accounting for one of the principal effects of a surface regression 
which varies along the length of the grain. This is the effect on tailoff 
characteristics which, in somewhat simplified terms, is treated in the 
present analysis by reducing the burning surface at tailoff in proportion to 
the length of propellant that has experienced burnout according to the linear 
model of the regressing surface. 

A discussion of the program input and output variables is presented in 
Section II of this report. This will serve as a guide to the user in proper 
specification of his problem and interpretation of the output. The discussion 
should also be read by those desiring to understand the mathematical analysis 
which is outlined as a step-by-step procedure in Section III. However, 

Section III may be omitted without loss of continuity by the user interested 
only in understanding how to specify his problem. Section IV presents the 
information required to operate the program and computer format and flow- 
charts. Once the input variables discussed in Section II have been specified, 
the information given in Section IV should suffice for anyone f ami liar with 
use of a computer having FORTRAN IV capability to obtain the program output . 
Also, the program listing and flowcharts in Section IV will be useful for 
determining how to modify the program and for a detailed understanding of 
the program logic. However, for those desiring to understand only the basic 
logic of the program, study of Section III is recommended, as it was written 
as a step-by-step procedure with the purpose in mind of conveying the basic 
program logic. Section V presents test cases and comparison with results of 
others, both theoretical and experimental. Again those seeking only to know 
how to use the program may omit this section, although the user will find, 
the example listings of input and output data, which are presented only here, 
helpful. Finally, the Appendix is a design data compilation resulting from 
use of the program which may be of aid to the designer in predicting effects 
of various design parameters on SRM characteristics and performance. 



II. DISCUSSION OF INPUT AND OUTPUT 


The computer program described in this report is capable of calculating 
performance for a large variety of different grains and combinations . On 
the other hand, to attain the goal of a simplified program, it is necessary 
to impose substantial restrictions on the geometric end constraints of the 
burning grain surfaces. Thus the accuracy of results will depend to a large 
degree on judicious selection of the geometric input variables. The difficul- 
ties encountered in approximating the geometric burning characteristics may 
be lessened, of course, by making use of the program options which permit 
specification of the burning geometry of part or all of the grain by using 
tabular inputs with the distance burned being the single argument. In any 
event, it is necessary that the program user have a clear understanding of 
the definition of each input variable so that he can appropriately specify each 
configuration to be analyzed. In addition, the user must understand how to 
specify proper design criteria and the physical characteristics of propellant 
and inert materials . 

In this section each input variable is defined in the order in which it 
is encountered in the program. The English or Greek symbol is given first 
followed by the computer symbol in parenthesis. Sketches of geometric 
characteristics are given for clarification. Where appropriate, additional 
discussion of the variable and recommended or typical numerical values are 
given. In addition, the present outputs of the program are defined. It is 
the aim of this section to provide a guide to the user in the specification 
of his problem and interpretation of the outputs . 


Propellant Characteristics 

Density of the solid propellant (slugs/in. 3 ). 


p p (HHO) 
a (A) 


n (H) 

ct eb (ALPHA) 


Propellant burning rate coefficient in the equation r = aP 11 
(in/sec-psia ). 

Burning rate exponent (l). 

Erosive burning coefficient in the Robillard-Lenoir burning 
rate rule (Equation III-ll)(in 2 * 8 -ft° ^/sec 1 ,e lbf° * 8 ) . 


0 (BETA) 
y (MU) 


Erosive burning pressure coefficient in the Robillard-Lenoir 
rule (1). 

Absolute viscosity of the chamber gases (lbf-sec/in-ft) . 


C* (CSTAR) Characteristic exhaust velocity (ft /sec). This is the gas 
dynamic characteristi c ve locity which can be obtained from 
the expression C* = t/RT /r where R is the gas constant and 
T c the adiabatic flame temperature of the propellant. 
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Primary Basic 
L (L) 

t w (TAU) 
v . (VCI) 

Ci 

V cf (VCF) 

D i * (DTI) 

a (ALFAK) 
n 

• e (theta) 

n 

V (BMP) 

*Ta (XT) 

z (ZO) 
o 


Motor Dimensions 


Total initial perforated grain length including gaps at slots 
(in.)* This is used only for calculation of the initial 
length Reynolds number and in the erosive burning rate equation. 

Web thickness of the controlling propellant length (in.)* (See 
ZO below for definition of controlling length) . If the grain 
is tapered, the length average value should be specified 
excluding lengths having additional taper not represented by 
ZO and segments located anywhere having a step_ decrease in thick- 
ness. Such step decreases must be handled by the additive 
tabular input option if they introduce two significantly 
different web thicknesses for the same grain type; e.g., two 
circular perforated grains . If a circular perforated grain 
is used, it is assumed that it will have the approximate 
average web thickness of the controlling propellant length. 

Initial volume of chamber gases (in 3 ). This is the volume of 
the empty case less the volume of the propellant . A reasonable 
estimate based on the case volume and the estimated volumetric 
loading density will suffice. 

Final volume of chamber gases (in 3 ). This is the volume of 
the empty case, A reasonable estimate will suffice. 

Initial diameter of the nozzle throat (in.). 

Exit half angle of the nozzle (radians). 

Cant angle of the nozzle with respect to the motor axis 
(radians) . 

Estimated length of grain at the aft end at start of first 
tailoff having an additional taper not represented by ZO or THETAG 
(in.). See Figure II-l. This variable permits the designer 
to specify an additional taper at the nozzle end of a circular 
perforated or star grain to produce regressivity shortly 
before tailoff. 

Difference in web thickness at ends of LTAP (in.). 

See Figure II-l. 

Initial difference between web thickness at the head and 

aft ends of controlling grain length, not including any 

initial length associated with LTAP or THETAG (in.). The 

controlling length of the grain is the axial length between 

the head end of the grain and the position of expected 

maximum Mach number in the port. It may be assumed, in most 

cases, that this length terminates whenever there is an abrupt 

decrease in web thickness near the aft end of the chamber or 

when 0 O is greater than 5°. 
u 




'Ph^ C ~ m + t0 fT f* men5 *° ns - ^Gci» ^Gni > and 0 G are used only for circular perforated groins 
The controlling grain length may he 1^ o? some lesser value depending on t^ posSi™ 

l°enX^l^ T rl\ r - iS tS ' le ”^ h —age -h thichness o^er tS g?^ UOn 
length excluding the lengths associated with I* a , e G and the head end dome. 
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Basic Performance 
Ay 2 (DELTAY) 

X out (X0UT > 

( AP/Ay ) t ( DPOUT ) 

C-, (2ETAP) 

Jt* 

(TB) 

\ (HB) 

X (GAM) 

E (RADER) 


Constants 


Initial desired burn increment (in.). This increment vill 
be used by the program for the initial 5 % of the web 
thickness burned, for the period shortly preceding and 
following tailoff and at such other times as the rate 
of pressure change is large. Larger increments will 
automatically be used at other times. An initial 
increment size of approximately 0.005t w is recommended 
but larger increments may be used for relatively flat 
burning surface versus distance burned traces. 

Distance burned at which the propellant breaks up (in.). 
This permits the option of specifying termination of 
burning resulting from possible structural breakup of 
the propellant. If this option is not desired, XOUT should 
be set to some large value; e.g., 1000 in. 

Rate of change of pressure with respect to distance burned 
at which the propellant is extinguished (lb/in 3 ). This 
permits the option of specifying termination of burning 
resulting from rapid depressurization. If this option is 
not desired, DPOUT should be set to some large value; 
e.g., 10,000 psia/in. However, if an abrupt tailoff is 

expected, a smaller value of DPOUT must be specified as 
the computer will not otherwise be able to handle the 
discontinuity. For large motors (120 in. dia. and up) a 
value of 500 psia/in is recommended and proportionately 
larger values for smaller motors. 

Thrust loss coefficient (l). In the absence of data to 
the contrary a value of 0.98 is recommended. 

Estimated burning time (sec). This and the variable that 
follows permits the program to calculate delivered specific 
impulse and thrust based on an assumed trajectory which was 
determined from analysis of typical large SRM applications. 

Estimated burnout altitude (ft). To obtain sea level 
performance characteristics, HB should be set equal to 
zero . 

Ratio of specific heats for propellant gases (l) , 

Radial erosion rate of the nozzle (in/sec). A value of 
0.0135 ips has been used by one propulsion company to 
predict the performance of 156 in. dia. motor with a 
graphite tape-phenolic impregnated throat . 
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Basic Properties for Weight Calculation 

(it ). (PIPK) Temperature sensitivity coefficient of pressure at constant 
pk K = A^/A* (1/°P). 

AT. (ETEMP) Maximum expected increase in initial grain temperature above 
design point conditions (°F). 

(ap/P) (SIGMAP) Statistical variation in the maximum calculated pressure(l) i 
i.e., one standard deviation in maximum pressure divided by 
the calculated value of maximum pressure. The value used is 
typically 0.03 for an untried design. Lover values may be 
used if supported by data on a proven design. 


o s /S (SIGMAS) 

yc 9 Statistical variation in yield strength of the case material 

(l). A value of 0.02 is recommended. Alternately Sy C may be 
taken as the minimum yield strength in which case SIGMAS is 0. 


n lS n 2 (Nl,N2) 


S 

ycnom 

(SYCN0M) 

D (DCC) 
cc 


Number of standard deviations in maximum chamber pressure and 
case yield strength, respectively, to be used as bases for 
design (l). A value of 3 gives statistical assurance that 
99 ‘ 5 % of the time the calculated limits of the variable will 
not be exceeded in tests. 

Nominal yield strength of the case material (lbs/in 2 ). See 
also SIGMAS. 

Estimated mean diameter of the motor case (in.). A rough 
estimation of the case and insulation thicknesses may be 
added to the grain diameter to obtain DCC. 


(PSIC) Safety factor on case thickness (l). This factor is super- 

imposed on the safety provided by use of the statistical 
limits on pressure and case strength. 


6 (DELC) 
c 


Specific weight of the case material (lbs/in 3 ). For steel 
the approximate value is 0.283 lbm/in 3 . 


L (LCC) 
cc 


Length of the cylindrical portion of the case including forward 
and aft segments (in.). 


n seg (HSBG) 


The number of case segments (l). NSEG is 1 for a mono- 
lithic case. 


h 


cn 


(HCN) 


Axial length of the nozzle closure (in.). See Figure II -1. 
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S (SYNNOM) 
yn 


Nominal yield strength, of the nozzle structural material 
(lb/in 2 ). 


(PSIS ,PSIA) 


K 1 ,K 2 (K1,K2) 


^ ins (PSIINS) 

6. (DELINS) 
ms 

K eh (KEH) 

K en (KEN) 

<s £ (dliner) 

T z (TAUL) 

W (WA) 
a 


Safety factor on the nozzle structural and ablative 
materials, respectively (l). No additional allowances are 
made in the computation for statistical variation, so PSIS 
and PSIA provide the only safety factors. Values of l.U 
and 2.0 are recommended for PSIS and PSIA, respectively. 

Empirical constants in the nozzle weight estimation equation 
(Eg. 63b). (lbm/in 3 , in/sec). For single, submerged conical 
nozzles using ablative exit cones with the steel body extend- 
ing over about one-third of the exit cone, the values recom- 
mended based on empirical analysis of large (120 to 156 in. 
dia. ) SRMs are K1 = 0,208 lbm/in 3 and K2 = 0.0925 in/sec. 
These values should also yield approximations of the weights 
for smaller and larger nozzles of the same type. For better 
results with larger or small motors or nozzles of different 
types, the K's should be redetermined. 

Safety factor on the nozzle insulation (l). No additional 
allowance is made in the computations for statistical varia- 
tion. A value of 2.0 has been used extensively in man -rated 
applications . 

Specific weight of the case insulation (lbm/in 3 ). DELINS is 
typically 0.0U62 lbm/in 3 . 

Erosion rate of the exposed insulation taken everywhere the 
same except at the nozzle closure (in/sec). KEH is typically 
0.003 in/ sec. 

Erosion rate of the exposed insulation in the nozzle closure . 
KEN is typically 0.012 in/sec. 

Specific weight of the liner (lbm/in 3 ). DLINER is typically 
0.033 lbm/in 3 . 

Thickness of the liner (in.). TAUL is typically 0.065 in. 

Any additional weight not considered elsewhere in weight 
calculations (ibm) . This includes weights associated with 
thrust termination and thrust reversal, thrust vectoring 
devices, thrust skirts (except for ”y” Joints), ignitors, 
instrumentation and final thermal protection of the motor 
case after burnout (see Eq. 6U , Section III). 



Input to Establish the Program and Basic Grain Configuration and Arrangement 


I (INPUT) 
op 


G (GRAIN) 
op 


S (STAR) 
op 


(NT) 

Q op (ORDER) 


1 for only tabular input . 

2 for only equation input . 

3 for a combination of 1 and 2. 

1 for an entirely circular perforated grain. 

2 for star grain only. 

3 for combination of 1 and 2. 

0 for an entirely circular perforated grain. 

1 for standard star (See Fig. II-2). 

2 for truncated star (See Fig. II-3). 

A standard star may not be combined with a truncated star. 

Number of thrust termination passageways in grain (l) 

NT is zero if there are no thrust termination passageways. 

1 if a star grain is at head end and a circular perforated 
grain at aft controlling end. 

2 if a circular perforated grain is at both ends. A star 
grain segment may still be present. 

3 if a circular perforated grain is at head end and a star 
at the aft controlling end. 

U if a star grain is at both ends. 

If grain = 1 , value of order must be 2 

If grain ~ 2, value of order must be 1+. 

It is important to realize that ORDER establishes the control- 
ling port area equations to be used. Thus if the nozzle end 
segment is not indeed the controlling one (the one that 
establishes the maximum Mach number in the port), ORDER should 
be specified to designate the actual controlling segment as 
the nozzle end segment. GRAIN, STAR, N1 and ORDER are not 
used for INPUT = 1, but values must be assigned for continuity 
of computer operations. 





gur -3. Truncated (slotted tube) star grain cross-section. The evolution of 
e urning perimeter is shown. Calculated variables are circled. 
First and third (final) burning zones are cross-hatched 



-13- 


Tabular Burning Surface and Port Areas (Not required for IHPUP = 2) 


\pT (ABPK > 
(fflSK! 

(AB1 ‘ K) 
V lAPHK> 


Burning area in the port ( in 2 ) . 

Burning area in the slots (in 2 ). 

Burning area at the nozzle end (in 2 ). 

Port area at the nozzle end of the controlling grain length 
(in 2 ). 


A pnT 


(APNK) 


Port area at the head end of the grain ( in 2 ) . 

APHK and APNK are not required when INPUT is 2 or 3; the 
equation inputs must be used to provide the information in 
these cases. Values of all A's used must be specified to 
completely describe the burning surface and port areas versus 
distance burned. The program computes intermediate values by 
linear interpolation. The number of values required is 
arbitrary and limited only by the storage capacity of the 
computer, but values must he specified for y = 0. Also, as 
an example of the procedure for specifying terminal values 
which must be followed, burning surfaces must be specified 
as zero at burnout of the tabular surfaces and at a 
y exceeding the highest anticipated calculated performance 
value. Separate input cards must be prepared for each value 
of y and arranged as described in Section IV of this report . 
The use of these tabular values in conjunction with the 
equation inputs (INPUT ~ 3) increases the flexibility of the 
program considerably. Frequently it is easy to estimate 
burning surface effects which end constraints on the equation 
inputs neglect. In this case a table of input values can 
be readily prepared from estimates of the effects. Also, the 
program can be run first without the tabular values and output 
used as an aid in obtaining the estimate. For example, the 
burning perimeters of a star grain can be determined in this 
way and the values used to estimate the effects of a head 
end closure on the star grain. 


Geometry for Circular Perforated Grain - Not required for completely star- 
shaped grain (GRAIN = 2) or for only tabular inputs (INPUT = l) 


D Q (DO) 
l> ± (DI) 


AD i (DELDI) 


Length average outside diameter of circular perforated 
grain, excluding lengths extending into the closures (in.). 

Length average inside initial diameter of circular perforated 
grain (in.). Only the controlling length excluding LTAP 
should he considered in the averaging 

Difference between the initial circular perforated grain 
diameter at the nozzle end of the controlling length and 
DI ( in . ) . 
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S (S) Number of burning flat ends of a circular perforated grain 

not including an extreme aft grain end (l). 

Angle burning surface element of circular perforated grain 
located at the extreme nozzle end of chamber makes with the 
motor axis (radians). See Figure II-l. THETAG must be set 
to zero if a star grain is located at the nozzle end (GRAIN = 3, 
ORDER = 3) or if aft end burning surfaces are represented by 
tabular values. THETAG is tt/2 if the circular perforated grain 
represented by equations has a flat burning surface located at 
the extreme nozzle end of the chamber. If THETAG is less than 
or equal to 5° , a value must be assigned (zero is satisfactory), 
but the effect of THETAG on burning surface area is not com- 
puted (See also LGNI and LTAP). THETAG is zero if the end 
surface is flat and inhibited. 

L Gci (LGCl) Initial total axial length of circular perforated grain 

represented by equation inputs not including gaps (in.). LGCI 
excludes lengths associated with THETAG. 

L Gni 


6 (THETCN) The contraction angle of a circular perforated grain bonded to 
' 1 the nozzle closure. See Figure II-l. Use an estimated value 

which yields approximately the correct volume of propellant 
burned. If a star shaped grain is located at the extreme nozzle 
end of the chamber or if tabular values are used to represent 
downstream burning surfaces, THETCN is tt/2. THETCN is also ir/2 
if the extreme aft end of the grain is inhibited, but only a 
flat ended, inhibited grain (THETAG = 0) which does not extend 
into the nozzle closure may be accurately represented. THETCN is 
zero for a burning flat end (THETAG = ir/2 ) which does not extend 
into the closure. 

The contraction angle of a circular perforated grain bonded to 
the head end. See Figure II-l. Use an estimated value which 
yields approximately the correct volume of propellant burned. 
THETCH is tt/2 if the extreme forward end of the circular perfor- 
ated grain (bonded or not) represented by equations is flat. A 
head, end flat burning surface is treated by proper specification 
of S. 


0 . (THETCH) 
ch 


Initial slant length of a burning conical circular perforated 
grain at the nozzle end (in.). LGNI is set equal to zero if 
THETAG is less than or equal to 5°. In this case the length 
otherwise associated with LGNI should be added to LGCI. If 
the error in burning surface area thus introduced is deemed 
significant, a correction may be introduced by making use of 
tabular inputs in combination with the equation inputs, Basic 
effect of small THETAG on tailoff may be accounted for by 
specification of LTAP. If a nozzle end burning surface is flat 
(THETAG = ir/2) , LGNI equals one half the difference between 
inside and outside local grain diameters. 


0 _ (THETAG) 

VJ 
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Basic Geometry for Star Grains - Not required for entirely circular perforated 
grains (GRAIN = l) or for only tabular inputs (INPUT - l) 


n (NS) 
s 


Number of burning flat end surfaces of a star grain not 
located at extreme nozzle end of the chamber (l). 


L Gsi (LGSI) 


n (NP) 
P 

R c (RC) 
f (FILL) 
n (NN) 


Initial total axial length of star-shaped perforated grain 
represented by equations (in.)* No provision comparable 
to the use of LGNI for circular perforated grains is made 
here to treat effects resulting from THETAG greater than 5°. 
Adjustments may be made, however, by use of tabular input 
values in conjunction with the equation inputs. Also, effects 
of taper, including additional small taper at the nozzle 
end, on tailoff may be treated by use of the variables ZO, XT 
and LTAP . 

The number of star points (l) . 

The star grain outside radius (in.). 

The fillet radius at star valleys (in.). See Figure 11-2,3. 

Number of burning flat end surfaces (0 or l) of a star grain 
located at the extreme nozzle end of the chamber (l). 


Special Geometry for Truncated Star - Not required for standard star grains 
(STAR = 1) 

R (RP) The length average initial radius of the truncation (in.). 

^ See Figure II-3. 

Tg (TAUS) The length average initial thickness of the propellant web 

at the bottom of the slots (in.). 


Special Geometry for Standard Star Grain - Not required for truncated star 
grains (STAR =2) 

T vg (TAUWS) The web thickness of a standard star grain (in.). See 

Figure II-2 . If the grain is tapered, the length average 
value should be used. TAUWS will be the same as TAU if 
the entire controlling length of grain is a star grain. 
Otherwise TAUWS may or may not be different from TAU. 

9^ (THETAF) Angular location of the fillet center of standard star 
from line of symmetry (radians) . 


6 (THETAP) 
P 


The apex angle of the star point (radians). 



Geometry of Thrust Termination Passageways - Not required if NT = 0. 


* TP (LTP) 


App (DTP) 

0 Tp (THETTP) 


Initial length of the termination passageway between the 
centers of gravity of perimeters of the bases (in.)- See 
Figure II-l . 

Initial diameter of the termination passage (in.). 

The acute angle between the axis of the passage and the 
motor axis (radians). 


T Teff (TAUEFF) Estimated effective web thickness at the termination port 

(in.). The user must Judge the distance burned at which the 
effect of the termination passage on modification of the 
burning surface gecanetry ceases to be significant. In general, 
this should be between two-thirds and full web thickness. The 
equation used to account for the burning surface is based on 
a passageway in a circular perforated grain terminated at 
the case by a flat inclined plane. Thus only a rough estimate 
of the effect of the termination passage is provided. 

Special Equation Inputs - Required only at the option of the user. May be 
used when INPUT = 1, 2, or 3. 


B^ (BBP) Additive burning surface input as function of Y for port 

burning surface (in 2 ). 

B^g (BBS) Additive burning surface input as function of Y for slot 

burning surface (in 2 ). 

B bn ( BBN ) Additive burning surface input as function of Y for nozzle 

end burning surface (in 2 ). In order to make use of the option 
of specifying the B's, a minor program modification is required. 
The B's are all set equal to zero in the present program. If 
this option is to be used the program statements assigning 
values to the B's are easily replaced with the desired 
equation inputs. 


Program Outputs 

The variables whose values are printed by the present program are 
defined below . Additional variables may also be printed with minor program 
modification. 



Outputs as functions of time 


t (T) 

y (Y) 

r n (RNOZ) 

r (KHKAD) 

P (PONOS) 
on 

P. (PHEAD) 
h 

l/J (PTAE) 

Mu (MJOZ) 

S G (SG) 

P a (PATM) 

A, + A. + A. 

Dp DS D 

(SUMAB) 

( CFVAC ) 

Fvac 

F (FVAC) 
vac 

F (F) 

Weight 

W _ (WP1) 

Pi 

W p2 (WP2) 


Operating time (sec.). This is calculated from the time 
initial equilibrium operating pressure is attained. 

Average distance burned (in.). 

Nozzle end controlling burning rate (in/sec). 

Head end burning rate (in/sec). 

Stagnation pressure at the nozzle end of the chamber (lb/in 2 ). 
Pressure at the head end of the chamber (lb/in 2 ). 

Port-to-throat cross-sectional area ratio (l) . 

Mach number at the nozzle end of the controlling grain length (l). 

Average perimeter of a star-shaped grain (in.). This will 
be printed as zero if there is no star grain or if tabular 
input values are used alone. 

The atmospheric pressure at altitude (lb/in 2 ). This is based 
on the assumed trajectory and specification of HB . 

The total burning surface of the propellant (in 2 ). 


The vacuum thrust coefficient (l). No adjustments are made 
in the calculations for nozzle divergence, cant, or other losses. 

The vacuum thrust (lbf). Adjustments are made in FVAC and F 
for nozzle divergence and cant, and other losses. 

The delivered thrust based on the assumed trajectory (lbf). 

Losses are included. 

and related values 


Propellant weight calculated from mass discharge rates (lbm) . 

Propellant weight calculated from the products of burning surfaces 
and incremental distances burned (lbm). 
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ma^ PHMAX ^ Maximum head end chamber pressure calculated by the program 
(lb/ in 2 ). 

W (WP) Arithmetic average of WP1 and WP2 (lbm). A check on the 

calculation accuracy is provided by comparison of this with 
WP1 or WP2. 

I S p (ISP) Average delivered specific impulse (lbf-sec/lbm) . Losses 

p are included as for FVAC and F. 

1^, (ITOT) Total delivered specific impulse (lbf-sec). Losses are 

included. 

I ( ISP VAC) Average vacuum specific impulse (lbf-sec/lbm). Losses are 

spvac included. 

P meop Maximum expected operating pressure (lb/in 2 ). This includes 

allowances for statistical variations and initial grain 
temperature . 

x cc (TAUCC) Thickness of the cylindrical section of the case. 

W c (WC) Total weight of the empty case (lbm). 

W n (WN) Total weight of the nozzle (lbm). 

W. (WINS) Total weight of the case insulation (lbm). 

ins 

W (WM) Total weight of the loaded case including added weight, 

WA (lbm). 


(ZETAM) Motor mass fraction; i.e., the ratio WP/WM (l). 

*T^M Total delivered impulse-to-motor weight ratio (lbf-sec/lbm) 



III. ANALYSIS 


This section of the report presents the mathematical analysis of the SRM 
design problem in a step-by-step procedure. The aim is twofold: l) ,to 

identify the fundamental physical relationships used and 2) to establish the 
basic mathematical logic of the computer program. 

No attempt has been made to justify completely the relationships used. 
The theory applied involves, in the main, basic internal ballistics and 
trigonometry applied in a more or less conventional way. The internal 
ballistics analysis follows one -dimensional frictionless compressible flow 
theory with the flow taken as reversible in the nozzle. Equilibrium pressure 
is established by iteration to find the chamber pressure which balances the 
mass flow rate relationship 


% - “o - 4<» c V /at 

A number of excellent references are available which explain the fundamental 
relations in some detail; e.g.. Reference 6, Chapters 5 and 6. Although in 
the present analysis fewer simplifying assumptions are made in the application 
of the physical principles than in the example and in similar references, 
the reader familiar with the basic internal ballistics of SKMs should have no 
difficulty in following the procedure and recognizing generalizations of the 
elementary approaches. In the development of the weight relationships, 
explanation of the assumptions made are given as changes in basic design 
criteria, new design approaches, or advances in the state-of-the-art may dictate 
changes to the equations used. 

Internal Ballistics 

At this point it is assumed that the port areas at the head and nozzle ends 
and burning surfaces areas are known as functions of the distance burned. 

The procedure for calculating these quantities is given later in this section. 
The computer program uses a subroutine in calculating the areas which is 
called as required. 


0. Set y = 0 and t = 0 

1. Compute A* and in turn from: 

a. A* = (hA)(D.# + 2E t) 2 
i n 


b. M n = (XA*/A pn ){2[l + (Y-l)M n 2 /2]/(Y+l)} (Y+l)/2(Y ~ l} 
where X = (A^ + ^/(A^ + + A^) 
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2 . 


3. 


k. 


5. 


6 . 


7. 

8 . 


9. 


10 . 


11 . 


12 . 


13. 


Compute 

u n = YC*M n {[2 /(y+ 1)] (y+1)/(y_1) /[1 + (T-DM n 2 /2]} !5 
P n / p on = [1 + (Y-l)M n 2 /2]" Y ^ Y_l) 
r = y J5 [2/(y+i)] ( y+i)/2(y_i) 


J = A* /A 

pn 

Compute first estimate of p on » 

P ob - IV'V + V + A ta )/A*] 1/<1 - n) (l + r2jV2)" /(:L ' n) 

Compute mass flow discharge rate through nozzle 
™D '■ A * P «/ C * 


P = (P /P ) P 
n n on on 


P. - P + 2in_u /(A , + A 
h n T) n ph ; 


pn 


Compute head end burning rate 

r = aP n 
h h 


) 


Compute nozzle (aft) end burning rate 


r n * °*n + “eb <V'V°- 8/L ° -2 ^VpVV 0 

The Robillard-Lenoir erosive burning rate rule is used. Reference 7 
provides support for this selection, but the relationship may be 
changed if desired to express r in teams of any of the variables avail- 
able at this point. n 

Compute the mass flow generation sate 


*0 * <V 2)[(r h + r » )( % + V * 2a C V 

If 1.001 m^ > m^ > 0.999 m^, the solution is complete for y=0. 

Go to step 16 


lU . If not , set * m Q > and 
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15 • P on = C * V** 

Return to step 8 

16. Compute the initial length Reynolds number 

R eL = % 

17. Set t = 0. 

18. Solve for M in 

e 

M e = A*/A e {2[1 + (Y-l)M e 2/2]/(Y+l)} (Y+l)/2(Y "' l) 

19. Compute 

a. P e /P 0 = [1 + (Y-l)M e 2 /2]“ Y/(Y " l) 

b. A = tfD 2/4 

e e 

20. Compute extimated altitude (Reference 8, pp. 7-13) 

h = 

21. Compute atmospheric pressure (Reference 8) 

P a = 11+.696 exp( -0.1+ 3103 x 10 -tf h) 
d.2.. Compute the thrust coefficients 

a. C F = r {2Y[l-(P e /P on ) (Y ‘ 1)/Y ]/( Y -l)}’ s 

+ (\/A*)[(P e /P on ) - 

“• °Fva= - C F + i 4 e«*)(V P on> 

23. Compute thrust 

23. a. P = C F (cos 6 n ) P on A*{c r (l + cos a^/2 

+ (l - cos “n )(A e / A*H( p e /P on ) " V P on ]/2 ^ 

23 ‘ b ‘ F vac * C P (cos V P on A * { C FVac (l+cos a n )/2 
+ (1-cos « n )U e /A*)(P e /P on )/2} 

The above thrust equations have been arranged so that the nozzle divergenc 
loss factor (1 + cos a n )/2 modifies only the momentum term of the thrust 
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/] Vc 


2U . Compute the instantaneous specific impulses and total delivered inpulses 

2b. a. I Bp = F/i^go 

2k. b. I = F 

spvac vac 

2 *. .. vn j4tj 

*• d ' Vac = J Vc J 4t J 

where the bars indicate here and below the arithmetic average value over 
Ay except at y^O where the initial value is used. 


Repeat steps 1 through 3 and 5 


25. Set y = A y^ and j = 1 

26. Compute a first estimate of new mass flow discharge rate using the value 

°? ^on beginning of the previous increment. 


■d * A * P on /C * 


Repeat steps 8 through 12 


27. a. Compute an estimate of dP^/dy, the rate of change of chamber pressure 
4P c /4y = (P aT/ ? av )(4i: av /4y) 

+ IP av^ A bp + V + V + V + V '<* 

where the subscript av indicates the arithmetic averages of nozzle and 
head end values. 

27. b. Compute the gaseous volume of the case 

V c “ V ci ♦ P A /j <% ♦ ♦ V>j 

28. If 1.002 ip > i 0 -jgfejsj [V. 4P,/Ay + + App ♦ A^)] > 0.998 ip. 


the solution is complete for y = £ Ay. 

j j 

29. If not, set 

“b * “<j - <Vay /12r2 C * 2) 4 V 4y > 

P onj = V W 


Go to step 31 


30. 


Repeat steps 8 through 12 and then 27 through 28 
31. a. Compute the time increment 


4t j “ 24y j /(? h + ? n> 
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31. b. Compute the time 

t = I at 

3 J 

32. Compute the difference between web thicknesses at the head and aft 
ends of the controlling grain length exclusive of any length near the 
nozzle end having an additional taper for the purpose of controlling 

tailoff characteristics (See discussion of t , z , L m and x„ , Section 
III). v o Ta Ta 


33. 


z 


Z At 

3 




3^ . Set j = 2 


35* If y >_ 0,05 t w and 


Repeat steps 18 through 2l(c. 


7 + Ay l < (t w " l z/2 i " *rJ 

,j-l 

l V c ? av /r2 C * 2) AP c /Ay lj-l = V 100 ’ 

set Ay^ = 10 Ay 1 

36. If not, set Ay = Ay 

J J- 

37 . Not used 

38. y = Z Ay 

J 

• a • Repeat steps 1 through 3 and then 2 6 through 33 and 33 through 38 
using A pn and A* in 1 

39. b. Set 3 = 3,*+,... Repeat steps 35 through 39b for each ,j 

until y >_ t w - 1 z/2 1 

39. c. Set z = z w and Ay = Ay^ (j is now w) 
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39. d. Compute the loss in length of turning propellant associated 
with additional tailoff controlling taper near the nozzle end of the 
grain . 

(!) lf = 0. set T £ = 0 

(2) If x_ > 0, set 

T * = - T » + V + 2/2 


If < 0, set = 0 

If T. > L m , set T„ = L_ 
l — Ta’ Jt T?a 


40. 

4o. 


If z <0 and 
w — 


a ‘ If l z w l > ^ y < t w + i z^/2 J , set 

% + V * V * Vo [1 + J 4j, / 2 w - iJr v /2z w ] 
-v 1 , 2 Ay j /z w - A V 2 V 


40. t. If |zj < Ay, or if y > T y + |z w /2|,set 

■^bp + 'Sjs + ^bn ~ 

41. If 7, > 0 and 

*V Z » + 4 V 2 2 v ] 
V 21 . 1 

^bp + ^Sjs + Vi " Vo 


4l. a. 


If z v J ^ and y < t y + z w /2, set 

V + + V = V [1 " 1 


0=v 


j 


+ \.+ ~ [ £ Ay./z - 

^>to j=v *y w 

4l. b. If z <_ Ay-, or y > t + z /2, set 
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1+2. 

a. 

Compute a first estimate of P during tailoff 

on 



P on ' t«P p + \ s + \„)M*] 1/(l - n) 

1+2. 

b. 

Set Ph = P 

on 

1+3. 

a. 

Compute - A*P q ^/C* 

1+3. 

b. 

4P c /iy * + \ s * \„>1 ♦ *6, + V )/iy 

1*3. 

c. 

r = aP n 
h on 

1*3. 

d. 

r n = r h 

1*1+ . 


m G r av + "^bs + ^bn^ 

1*5- 

Not 

used 

1*6. 


Repeat steps 27b through 29, then 47a, b, and c below as required 



until 28 is satisfied 

1*7. 

a. 

P . * P , . _ v + ( AP /Ay) Ay, 
onj on(j-l) c 

1*7- 

b * 

P = P 
h on 

1*7. 

c. 

Repeat steps 1+ 3b through 1+6 

1*8. 

a. 

Repeat steps 31a through 33, 36 and 38 

1*8. 

b. 

Repeat steps 40 through l+7c for other values of j > w until 



P <30 psia, or 
on — * 

|AP c /Ay| > |(AP/Ay) out |, or 



y - * out ’ or 



( % + A bs + V = ° 


If 

any of the conditions in step 48b are met , begin "half-second trace 

1*9. 

a. 

Set p ht = \ 

1+9. 

b. 

Set t. = t 
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*»9. c. P n - P ht exp[-r 2 C* 2 12(t - t t )/V c ] 

49. d. Set P = P 

on n 

1+9- e. Compute = P A*/C* 

49. f. Compute t - Z At. 

J J 

49. g. Set At ' = 0.5 sec. 

49. h. Repeat steps 49c through. 49g ten times or until 

1 . 3° P sia 

50. Not used 

51. a. Calculate weight of propellant "burned 




W . = g I 
px *0 J 


51. 

b. 

^p2 P p®o 

] ( v + 

51. 

c. 

W = (w + 
P Pi 

OJ 

Cvj 

Cl 

* 

52. 

a. 

Compute 




V = X T /W P 

52. 

b. 

I = 

spav "vac 

I- /w 

T vac p 


'J 


53. Determine the maximum head end pressure P, 


max 


Inert Weight Calculations 

54. Calculate the maximum expected instantaneous pressure at the upper 
limit of initial grain temperature 


meop “ P h max [l + n l (o P /P)] ex P [(n p V T i ] 
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55. Calculate the minimum expected case material yield strength 

S - S [1 - n„ (o_ /S )] 

yc yc nom 2 S yc 

yc 

56. Compute the nominal thickness of the cylindrical portion of the case 

T lit P D /2S 

cc c meop cc yc 

57 • Compute the weight of the case material in the cylindrical length 
including forward and aft segment joints 

W — 7tt D <5 L [1+ (n - l)(ltOx /L )] 

cc cc cc c cc seg cc cc 

Note: A joint is considered ^ nominal case thicknesses thick and 

10 nominal case thicknesses long. 

58. Compute the hoop-stress thickness requirement for the head end and aft 
end closures. 

t = lj; P D /ItS 
l c a r c meop cc yc 

59. Compute the weight of the head end dome 

W . = 2.5 ttD 2 6 t ,/2 
ch cc c cd 

Note: The 2.5 multiplying factor is an allowance for the igniter boss 

and forward thrust skirt attachment joint ("Y" joint). A similar but 
larger (U . 5 ) allowance is made in the aft closure weight calculation for 
the nozzle and aft thrust skirt attachment joint. In case only one skirt 
is used, the dual allowance should still roughly approximate the net 
weight requirement because of the greater load carrying requirement placed 
on the single thrust reacting element. 

60. Compute the estimated weight of the nozzle closure 

W » U.5 irD h x , 6 /2 
cn cc cn cd c 


6 l. Compute the total estimated weight of the motor case 


W = W +W.+W 
c cc ch cn 
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62. Compute the nozzle expansion ratio 

e = A /A* 
e 

63 . Compute the estimated weight of the nozzle 

W n = (K 2 t f * a + (c-A P meop D* * 3 / 3 ^ 1/11 + (sin a n )/2] 

64. Compute the estimated weight of the case insulation 

W. = t_ \p. <5. ttD {K [0.80 D/2 

ins f r ins ms cc 1 en 1 cc' 

+ (s+n ) t /2 
s w' 


+ 0.15 (L - t S - t n )/i|). ] 

cc w w s ins 

+ K 0.80 h } 
en cn J 


Note: Head and nozzle end domes are assumed exposed for 80 percent of t^> 

The cylindrical length, exclusive of the areas exposed by circumferential 
slot burning, is assumed exposed for 15 percent of t-^ with no safety 
factor. No allowance is made for final thermal protection of exposed 
areas for which an allowance may be included in W a . 


65 . Compute the estimated weight of the liner 

W = t. 6 ttD [D J2 + L + h ] 
t & 1 cc 1 cc' cc cn 

66. Compute the estimated total motor weight 


W = W + W +W. +W.+W + W 

men ins tap 

6 7- Compute the motor mass fraction 


! M - V W M 

68. Compute the total impulse to motor weight ratio I^/Wy 
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Burning Surfaces and Port Areas 


The geometric relations described below are treated in the computer 
program as a major subroutine. It should be noted that end effects are 
treated by approximations both in the specification of input variables and 
in several instances in the mathematical relationships used. 


AO. a. 


AO.b . 


AO .c . 


Select an option for the type of input. 

If only tabular input , set I = 1 

If only equation input , set I = 2 

If a combination input , set I = 3 

Select an option for the grain configuration. 

If only a circular perforated grain, set G = 1 

If only a star grain, set G =2 

op 

If a combination grain, set G =3 

op 

Select an option for the type of star. 


If a standard star, set S =1 

op 

If a truncated star, set S =2 

op 

AO.d. Select an option for the grain arrangement. 

If a star grain is at head end of the chamber and a 

circular perforated grain at nozzle end, set Q =1 

op 

If circular perforated grains are at both ends of 

the chamber, set Q =2 

op 

If a circular perforated grain is at the head end of the 

chamber and a star grain is at the nozzle end, set Q =3 

op 

If star grains are at both ends of the chamber, set Q = U 

* op 

AO.e. Specify the number (rO of thrust termination passageways in the 
grain. 


AO.f. 


Specify tabular input values of burning surface areas 
and A^ nT versus y. Set each equal to zero if I = 2. 


^bpT’ ^bsT 


> 
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AO.g. 

Al. 

A2 , 
A3, a. 
A3.b. 

Ah. 

A^ . a . 

AU.b. 


Ak.c. 


If I - 1, specify the head end and nozzle end controlling 

port°§reas A m versus y. 
pnT 


If I op = 1, set AA bT , 


If I = 2 or 3 

op 

If G = 1 or 3 
op 

If G = 2, set A. 
op * Dpc 


^bpc %s \ sc ^bss ^bnc ~ "^b ns ® 

Go to step A51 

Go to step A3a 

Go to step AU 

■*bsc ~ ^bnc - ® 


Go to step AT 

Compute the burning surface areas of tubular grain segments: 

For slots 

Sc ' ’S/[B o 2 -(B 1 + 2y) 2 ]A for \ ec > 0 
A, =0 for A, <0 

DSC DSC — 

For perforations 

If 2y + Ih > D q} set A (jnc = A^j = 0 Go to step A 5 

If not and if 0 Q £ 5° (See Figure II-l) compute 

L Gc = L Ccl - y (oc,t S 4 oot V 81,4 

Sc “ * tD i + - T t - Syl for Sc ’ 0 

A, =0 for A, < 0, or 
Dpc Dpc — ’ 

If 0 Q > 5°, compute 


Sc = ’ (B i 4 2y)lL 0cl - y cot 9 ch - T * 

- [s + tan ,(0G/2)]y} for > 0 

A. =0 for A, < 0 
Dpc Dpc — 

For the nozzle end surface 

If 0_, < 5°, set A. =0 
G — * Dnc 
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If 6„ > 5°, compute 

U 

\nc " [L Oni - y cot (9 G + e on> 

- y tan (0 g / 2)][D ± + + y 

+ L„ . sin 0^ + y esc (0^ + 9 ) sin 0 ] for A, > 0 

Gni G ^ G cn cn Hanc 

Note: the above equation represents the entire surface as conical 

although a non-conical surface area evolves if 0_ + 9 > tt/2. 

& G cn 

A. =0 for A, <0 
Dnc one — 

Compute the controlling port areas 


If If Q = 1 or k 
op 

If Q = 2 or 3, compute 
op 

A , = tt[D. +2 I r. . At.] 2 /!* for A , < ttD 2 /l* 

ph x j hj j J ph — o 

A , = ttD 2 A for A , > nD 2 /l* 
ph o ph o 

If Q = 3 or k 
op 

If Q = 1 or 2, compute 


Go to step A5b. 


Go to step A 6 


A = u[D. + AD. +21 r . At.] 2 /!* for A < irD 2 /U 
pn l i j nj y pn — o 


A = irD 2 /k for A > irD 2 /U 
pn o pn o' 


a6. 

If G 

= 1, 


op 


A7*a. 

If G 

= 2 or 3 and S 

op 

op 

A7.b. 

If G 

= 2 or 3 and S 

op 

op 

A8. 

Compute 

the star point 1 


e 

= Tt/n 


s 

P 

A9. 

Compute 

the approximate 


a 

s 

= E - x - f - : 
c s 


Go to step A50 
Go to step A2Q 
Go to step A8 



-32- 


A10. 

All. 
All. a. 

All.b. 

A12. 

A13.a. 

A13.b. 

Aik. 

A15- 
Ai6. 
AIT ,18 


Compute the angular location of the "burning slot sides 

G = 0. - arc sin [ (f+y)/(R +y) 3 
Si s p 

Compute the perimeter of the grain. (if 0 g ^ < 0, Go to step A22 ) 

If y <_ r , compute 
" s 

S 0 ■ 2n p U K p + y) 9 sl + ‘s - t(R p + y) cos (9 s • 9 sl ) 

- R p 3 + ir(f+y)/2} Go to step Al^ 


If not , compute 


Q n = arc sin {[B o 2 -U c + Rj 2 -(f+y) 2 ]/2(f+y)U s + B ) } 


C ' 1 ' C S p 

If 6 0, compute 


S G = 2 ”p' (B p + y) 9 si + *s - [(R p + y) 008 {6 s ’ 9 sl ) 


- R p ] + (y+f) e c } 

If 0 <0 and y < B -R , compute 

c c p 

S G ’ 2n p 6 sl + l V - (f+y)21 ' 

- t(a +y) 2 - (f+y) 2 ]* 5 } 


Go to step Aik 


Go to step Aik 


If 0 <0 and y > E -R , set S- = 0 
c * — c p’ G 

Compute the initial "burning area of flat ends of grain. 


A Gs = ' (B o 2 -V ) ■ “p ( ’ f2/2 + 2 t s f> 

Compute burning slot end surface areas 
A bss = n s ^ A Gs ‘ J ^Gj r avj 


Set 




ss 


= 0 if 


< o or S G = 0 


Compute burning flat nozzle end surface area 

A, = n [A n - £ § r . At ] 

Dns n Gs j G,} nj j 


Set A. „ 
ons 

>19 Not used. 


= 0 if 


A, '< 0 or S. = 0 
ons G 


Go to step A31 



-33- 


A20 . 


A21 . 


A22 . 
A22 .a. 


A22.t>. 


A22 .c . 


A22.d. 


A23. 


Compute the half angle of the standard star grain (See Fig. II-2). 

0 = */n 

s p 

Compute the radial position of the fillet center 

£ = R - t - f 

few 

Compute the grain perimeter (Equation A22a and A22b are taken from 
Reference 6, p. 303) 

If (y+f)/£ < (sin 0 J/cos (9 /2) and y < t , 

I I p V 

S G = 2n p £ f {(sin9 f )/sin (0 p /2) + [(y+f )/£ f 3 [ff/2 + 0 g 

- 0/2 cot (0/2)] + 0 - 0 } Go to step A23 

ir Sr & X 

If (y+f)/t f <_ (sin0 f )/cos(0 p /2) and y >_ x^, 

[R c 2 -£ /-(y+f)*] 

S = 2n {(y+f) arc sin{ - (0 C - 0/2)} 

0 p 2* f (y + f) 6 p 

+ U sin9 )/sin(0 /2) - (y+f) cot (0/2)} 
tip P 

If (y+f)/£ > ( sin0 )/ cos ( 6/2) and y < x , 

lip W 

S G “ 2n p £ f l^y +f )/ £ f He s + arc sin (— sin 0 f ) ] + 9 g - 9 f } 

Go to step A23 

If (y+f)/*- > (sin0 )/cos (0 /2 ) and y > x , 

lip — V 

r l t 

S G = 2n,p (y+f) {e f + arc sin ( f — sin 0/ 

[R c 2_£ f 2 " ^ +f)2] 

~ arc cos } for S > 0 

2i f (y+f) 

S„ = 0 for S„ < 0 
u u — 

Compute the initial burning areas of slots and flat nozzle end. 

(If 6 ^ < 0, Go to step A1.1Q 

Aq s = itR c 2 - n// {sin 0 f (cos e f - (sin 0 f ) cot (0/2)] 

+ 0-0 + (2f / l ) [(sin 0 )/sin (0 /2) 

s f f , f P Q 

+ 0 g - 0 f + (f/2t f )(Ti/2 + e 5 - e p /2 - cot -|)]} 
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A24. Compute the burning slot surfaces 

\ss * n s [A 0s * \ % ? avj f ° r * 0 

J 


= 0 for 




A25. Compute the burning flat nozzle end surface 

A, = n [A-, - £ r , At,] for A, >0 

Dns n ‘•TJs j Gj nj y Dns 

A. = 0 for A. < 0 
ons ons — 

A26-30 Not used 

A31. Compute the port burning surface areas 

A, = [L . - T, - y (n +n )] S_ for A^ >0 
bps Gsi Z J s n G Dps 

A, =0 for A, <0 
Dps Dps — 

A32-35 Not used 

A36. Compute the controlling port areas of the grain segments 


A36.a. If Q = 2 or 3 
op 


Go to step A36b 


If Q = 1 or It, compute 

op 

A = ttR 2 - A_ + 
ph c Gs 


S„, r. . At. for A . < irR 2 
Gj hj j ph c 


A . = irR 1 for A , > ttR z or S_ < 0 
ph c ph — c G — 


A36.b. If Q = 1 or 2 
op 


Go to step A50 


If = 3 or 4, compute 

A = irR 2 - A„ + I 3_, r , At, for A < irR 2 
pn c Gs j Gj nj j pn c 


A = irR 2 for A > irR 2 or S_ < 0 
pn c pn — c G — 
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A37-U9 

A50.a. 

A50.b. 


A51. 

A52. 

A52.a. 


A53. 


Not used 


If there are no thrust termination passageways in the grain (N,j=0), 


If n T f 0, compute change in A^ 


Go to step A51 


AA bT = V [(l) Tp + 2y)(£ Tp' y/sin 6 Tp ) 

- (typp + 2y) 2 /** * < y' fD Tp , ' 2 )< D T p/ 2)(1 " 1/sin e Tp)l for y i T Teff 


iA^ = ° for y > x Teff 


Compute special equation surface area inputs , B^ g , and B^ n 


Note: the B's are all set equal to zero in the present program hut 

the zeros may he replaced by any functions of y desired simply by 
changing three program statements giving the B's 


Compute the total surface areas 

^bp ^bpT + "^bpc + ^bps + A-A bT + ^bp 

^bs ^bsT + ^bsc + "^bss + ^bs 

■^bn ^SanT + \>nc + ^bns + B bn 


Note: the present analysis treats A^p and A^, s identically. They 

are kept separate, however, to permit possible alternate treatment. 

Compute the burning surface areas applicable to tailoff surface 
area calculations 


If z >0, set 

w — 


*tm = (A bp + + V 


7 ~ 


i Ay /2 

J=w ° 


^to = + \s + 


y + 


If z. < 0, set 

W 


(z /2) - 

w 


J 

2 Ay ,/2 

J-w J 


■^bm " ^\p + ^s + ^bn^ j y 
^bto _ ^^bp + "Sds + ^bn^ |y 


(z/2) - E Ay J2 
i= w J 

J 

E Ay /2 
j=w J 



IV. THE COMPUTER PROGRAM 


This section contains the instructions for preparation and arrangement 
of the data cards. Also a complete listing of the program statement is 
given followed by flowcharts showing the sequence of computer operations . 


Data Card Usage 

The data formats have been established to allow the operator to look 
at the card and know which variables are represented on it. The format 
should be followed to insure correct reading of the inputs » The various data 
cards are as follows: (See Table IV-1 for the complete format) 



I. 

Number of configurations 


II. 

Initial zero values 


III. 

Propellant data 

IV 

a&b. 

Motor geometry (2 cards) 

V 

a&b . 

Performance data (2 cards) 


VI. 

Input, grain, etc. 

VII 

a&b. 

Data for tabular input (2 cards) 

VIII 

a&b. 

Tubular grain input (2 cards) 


IX. 

General star data 


X. 

Truncated star data 


XI. 

Standard star data 


XII. 

Termination port data 

XIII 

a-e . 

Data for weight calculations (5 cards) 


Cards II “VI must accompany each, and every configuration even if only 
one parameter changes. If INPUT = 1 or 3, a number of sets of data cards 
VII a&b are used. Consider a test run on three different configurations: 


a. A combination grain with only tabular inputs. 

b. A combination circular perforated and standard star grain with 
only equation inputs (no termination ports). 

c. A truncated star grain with two termination ports using both 
tabular and equation inputs. 


- 36 - 
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The correct order for the data cards is: 


Configuration a 

Configuration b 

Configuration c 

i. 

I 

12. 

II 

21. 

II 

2. 

II 

13. 

III 

22. 

III 

3. 

III 

lit. 

IV a&b 

23. 

IV a&b 

b. 

IV a&b 

15. 

V a&b 

2b. 

V a&b 

5. 

V a&b 

1 6 . 

VI 

25. 

VI 

6. 

VI 

17. 

VIII a&b 

26 . 

VII a&b (y=0) 

7. 

VII a&b (y=0) 

18 . 

IX 

27. 

IX 

8. 

VII a&b (y=y 2 ) 

19. 

XI 

28. 

X 

9. 

VII a&b (y=y 2 ) 

20. 

XIII a-e 

29. 

XII 

10. 

11. 

VII a&b ( y=y 3 ) 
XIII a-e 



30. 

31. 

VII a&b (y= yi ) 
VII a&b (y=y 2 ) 





32. 

VII a&b (y=y 3 ) 





33. 

XIII a-e 


Note that when tabular inputs are used alone, as in Configuration a., the 
sets of cards VII a&b (h sets for Configuration a) . giving the tabular inputs 
follow each other sequentially. When combination tabular and equation inputs 
are used, as in Configuration c., the sets of cards VII a&b for y=0 follows 
card VI and the remaining sets of cards VII a&b follow card XII or the next 
lower applicable number if XII does not apply. Note also that if a particular 
datacard does not apply, no card is used. For example, X is omitted for 
Configuration b. and XI for Configuration c. 

Program Listing 

A complete program listing is presented in Table IV-2 
Flowcharts 


IV-3. 


Flowcharts of the computer program are presented in Figures IV-1 through 



TABLE IV-1 


DATA CARD FORMATS 


ccl cclO cc20 cc30 ccl+O cc50 cc60 cc70 cc80 


NUMBER OF CONFIGURATIONS TO BE TESTED =— 


i | I H 

0.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.0 


RHO— .- 


A=-. N— .- 


Ill 

ALPHA= — . - BETA= MU=+-.— — E+— CSTAR= 


IVa 

L= . TAU= — . — VCI=+- . E+ — VCF=+-. E+ — j DE= .— ;DTI= .— 


IVb 

THETA— . ALFAN— . — — LTAP= .— XI=— .— 20=—.— 


Va 


DELTAY— . - — Z0UT= .— DPOUT=- 


j RADER=-i 


ZETAF=. TB= .- HB= . GAM—. 


Vb 


VI 


INPUT=— GRAIN=— STAR= — NT= . ORDER=- 


VI la 

Y= ABP0RT=+- . E+ — ABSL0T=+- . -E+— ABN0Z=+-. E+- 


Vllb 

APHEAD=+- . -E+— APN0Z=+- . E+- 
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TABLE IV-1 (Cont'd) 


GCl CClO 


cc20 


CC30 


ccl*0 cc50 cc60 


cc70 


Villa 

D0= . DI=---. DELDI- . S= L THETAG=- . • 


I Vlllb i 

LGCI= — - — . — LGNI= THETCN-. THETCH=- .• 


NS=- 


! EP=- 


LGSI=- 


IX 

NP=— . BC= — - . FILLET= HN-- 


PIPK=.— 


I 


H2=— 


LCC=- 


TAUS= 


| THETAF=I. THETAP=-.- 


XI 

TAUWS “ — .• 


XII 


HTEMP=- 


Xllla 

[SIGMAP-. SIGMAS=-.- 


SYCNOM= 


HSEG= — . HCN=- 


Xlllb 

DCC= PSIC= — .■ 


XlIIc 

SYNWOM=- 


PSIA= — . — Kl=-. 


KEH=- . KEN=-.- 


XHId 
K2=- . PSIINS=- 


LTP= — - . — DTP= — . — THETTP=- . TATJEFP= — 


iSfl=- 


DELC=- 


f 

PSIS=- 


DELINS= 


XI He 

DLINEB=- . TAUL=- 


WA=~ 
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cc80 


! 



TABLE IV-2 


Gfc l 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


********************************************************************* 


* SR M DESIGN ANC PERFORMANCE ANALYSIS * 

* PREPARED AT AUBlRN UNIVERSITY 4 

4 UNDER MCD. NO. 6 TO COOPERATIVE AGREEMENT h I TH * 

* NASA MARSHALL SPACE FLIGHT CENTER * 

* * 

* ANALYSIS BY R.H. SFCRZINI * 

* PROGRAMMING BY J.M. LYON * 

* AEROSPACE ENGINEERING DEPARTMENT * 

* AUGUST 1972 * 


********************************************************************* 
REAL MGEN,KCIS,MNOZ,MM,JRCCK,N,L»MEl»ME,ISP, I TOT , MU, MASS, ISP VAC 
REAL M,N2,NS£G v Kl»K2,KEH,KEN,NS,LCC,LTAP 
REAL M2, MCBAR, ISP2, ITVAC 
COMMON/CCNSTl/Zfc,AE,AT,THETA, ALFAN,G 
CCMMCN/CCNST2/CAPGAM,Mt , BCT ,Z E T AF , TB , HB , GAM 
CCMMCN/CCNST3/S,NS 

CCMMCN/VAR I Al/Y, T , CELY , CE L T AT , PCNCZ » PHE AO , RNCZ , R HE AD , SUMAB , PHMAX 
COMMGN/VARl A 2/ A E PORT ,ABSLOT, ABNCZ, APHE AC , APNCZ , D AO Y, ABP2 , ABN 2 , ABS2 
COMMON /VAR I A3/ IT CT , ITVAC, JROCK , ISP* I SPVAC , MC I S , MNCZ , SG ,$UMMT 
CCMMCN/VAR I A4/ RN T , RHT , SLM2 ,R I, R2,R3,RHAVE ,RNAVE , RB AR , Y8 , KGUNT ,TL 
COMMCN/VAR1 A5/AEMA IN, ABTO, SUMCY 
RfcAC ( 5 , 5CO ) NRLNS 

**4**4*4*******444********************4****************************4* 

* READ IN THE NUMBER CF CC NF I GURA T I CN S TC BE TESTED * 

*************** 444***4*44* ******************************************* 

5GC FORMAT ( 4 2X * 12) 

CC 901 1=1, NRLNS 
VtR I TE ( 6 , 602 I I 

602 F0RMAT(1F1,42X, ‘CONFIGURATION NUMBER «,I2) 

RE AC (5, A 99) SUMC Y , ANS ♦ Z M , Y ,T , CELT AT , RNOZ , RHEAC , SUM A8, PHMAX , SUK2 , I T 
1CT,RHT ,R NT, Rl,R2 ,R3, RFAVE,RNAVE, REAR , ITVAC, SUMMT 
*<******************************************************************* 

* SET INITIAL VALUES CF SELECTED VARIABLES ECUAL TC ZERC 4 

* * 4 *N0 T E 4 ** THESE VALUES MUST BE ZERCEC AT THE BEGINNING CF 4 

* EACH CCNF IGLRATICN RUN * 

***************** ********************4$*******,* *** ******************* 

499 FCRMAU22F3.il 

REAC (5,501) RFC,A,N,ALPHA,B£TA,Mt,CSTAR 
***************** ******************* ********************************* 


4 REAC IN BASIC PROPELLANT CHARACTERISTICS * 
* * 
4 RHC IS THE DENSITY CF THE PROPELLANT IN SLUGS/IN4*3 4 
* A IS THE BURNING RATE COEFFICIENT * 
4 N IS THE BURNING RATE EXPONENT * 
4 ALPHA ANC BETA ARE THE CONSTANTS IN THE EROSIVE BURNING 4 
4 RELATION OF RCBILLARD AND LENIOR * 
4 MU IS THE VISCOSITY CF THE PROPELLANT GASES * 


/ 
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TABLE IV-2 (Cont'd) 


* CSTAR IS TEE CHARAC TERl ST IC EXHAUST VELOCITY IN FT/SEC ♦ 

501 FORK AT < 4 X , F8 . fc , 3 X, F 6 .4 , 3X , f 5 . 3 , 7X,F4.1,6X,F5. l»4X,E11.4 t 7X,F5.0) 

UR I TL (6,603) RFC,A,N,ALPHA,eETA, MU, CSTAR 

603 FORMAT ( / / , 20X , * PRC P EL L ANT CHARACTER 1ST ICS 13X , *RHC = »,F8.6,/,1 
1 3X , • A= ' ,F6.4,/, 13X, 'N= • ,F5.3, /, 13X, »ALPFA = » , F 4 . 1 , / ♦ 1 3 X , 1 BE T A = • 
2,F5.l»/, 1 3x » * M L- • ,1PE11.4,/,13X, »CSTAR= *,1PE11.4) 

RE AD (5, 5 02) L , T A l , VC I , VC F , CE, CT 1 , THE TA , A LFAN » L TAP , XT » ZC 

*$************<******* 4 ********$itf***1f:**m:*+1+ + * + 4* + 

* READ IN BASIC MOTOR DIMENSIONS * 

* * 

* L IS THE- TCTAL LENGTH OF THE GRAIN IN INCHES * 

* TAU IS THE AVERAGE WEB THICKNESS OF THE CONTROLLING GRAIN * 

* LENGTH IN INCHES * 

* VCI IS THE INITIAL CHAMBER VOLUME IN IN**3 * 

* VCF IS THE FINAL CHAMBER VOLUME IN IN**3 * 

* CE IS THE CIA METER OF THE NOZZLE EXIT IN INCHES * 

* DTI IS THE INITIAL DIAMETER CF THE NCZZLE THROAT IN INCHES * 

* THETA IS THE CANT ANGLE CF THE NCZZLE WITH RESPECT TC THE * 

* MOTOR AXIS IN RADIANS * 

* ALFAN IS THE EXIT HALF ANGLE CF THE NOZZLE IN RADIANS * 

* LTAP IS THE LENGTH CF THE GRAIN AT THE NCZZLE END HAVING * 

* ADDITIONAL TAPER NCT REPRESENTED BY ZC IN INCHES * 

* XT IS THE DIFFERENCE IN WEB THICKNESS ASSOCIATED WITH LTAP * 

* ZO IS THE INITIAL DIFFERENCE BETWEEN WEB THICKNESSES AT THE * 

* HE AC AND AFT ENDS CF THE CONTROLLING GRAIN LENGTH * 

$$¥$***#***$**** tf*****^***^*****#*;*#*****#********, $*,},***$$ *$****4,*.. 

502 FORMAT (2X,F5.0t5X,F5.2,5X,EU.4,5X,EH.4,4X*F6.2,5X f F6.2,/ t l4X»F6. 
14,7X,F6.4,7X,F6.2,5X,F5.2,5X,F5.2> 

WRITE ( fc,604 ) L » T AU ,VC I , VCF , C E , D T I , THE T A f AL F AN , LT AP , XT , ZC 

604 FORMAT (/ / ,2CX, 'BASIC MOTOR DIMENSIONS' »/,13X,*L= • , F 5 « C , / , 13 X , » T AL 
1= • ,F5.2,/, 13X, • VCI= * ♦ IPE 1 1.4,/, I3X, 'VCF= • * IPE 11.4,/,I3X,'0E= 
21PE11.4,/,13X, »CTI = ' , IPEI1.4,/, I3X, ' TH6TA= ' , IPE i 1 .4, / , 13X , • ALFAN = 
3 • , 1PEU .4,/, I3X,'LTAP= • , IPE 11 .4, /, 13X f * XT* « , 1 PE 1 1 . 4 , / , 1 3X , * ZG= 
4* T IPE 11. 4 ) 

READ! 5,503 ) DELTA Y,XGLT ,0 PCUT , ZET A F ,TB , H8 , GAM , RADER 
*44******«*****4***4*************************4***********4**4******** 

* READ IN BASIC PERFORMANCE CONSTANTS * 


CELTAY IS THE DESIRED BURN INCREMENT DURING TAILOFF IN INCHES * 
XOUT IS THE DISTANCE BURNED IN INCHES AT WHICH THE PROPELLANT * 
BREAKS UP 4 

DPCUT IS THE DEPRESSURIZATION RATE IN Le/IN**3 AT WHICH THE * 
PROPELLANT IS EXTINGUISHED * 

ZETAF IS THE THRUST LOSS COEFFICIENT 4 

TB IS THE ESTIMATED BURN TIME IN SECONDS * 

HB IS THE ESTIMATED BURNOUT ALTITUDE IN FEET * 

GAM IS THE RATIO OF SPECIFIC HEATS FCR THE PROPELLANT GASES * 



AGE 
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TABLE IV- 2 (Cont'd) 


C * RADER IS THE RACIAL EROSION RATE OF THE NOZZLE THROAT IN * 

C * INCHCS/SEC * 

C ’M***44*'M*44*4**444'****4»*******4************>e‘*******4**44******4t** 

503 FGRMAT (7X,f 5.3,6X,F7.2, 7X,F7.2,7X,FA.3,AX,F5. l,AX,F7.C,5X,F5.3,/,e 
1 X » F6 ♦ A ) 

WR I TE ( 6» 606 ) CEL TAY , XCl T , DPCUT , Z E TAF » TB , HB t GAM , RADER 
606 FORMAT (//, 15X, 'BAS IC PERFORMANCE CCNST ANT S • t / # 13X , ' DEL TAY= ' ,F5.3, 
I / » 13 X » ' XCLT= • ,F7.2,/,13X, ‘DPOUT* S F7 . 2 f / , 13 X , • ZETAF= *,F5.3,/,13 
2X,*TB= • ,F5.1,/ ,13X, »FB = • , F7 .0 , / , 13X , • GAK = 1 »F5.3»/»13X»*RA0ER= • 
3,F6.A) 

HNl=.85 
ME 1 = 7 . 0 
G = 32 . 2 
Z*ZO 
S = 0.0 
NS=0 . 0 
KCUN f = 0 
ABMA IN=0 *0 
AB TC=C .0 
DELY=CELTAY 
TCP=GAM+1. 

BOT=GAM- 1. 

VC=VCI 

ZAP=TCP/ (2.4BCT ) 

CAFGAM=SCRT (GAM )*( 2./TCP )**ZAP 
AE=3.1A159*CE**2/A. 

1 IF (XT. LE. 0.0) 7 L = Q .0 
IF (XT.LE.C.O) GC TO AC 
TL=(Y-TAU+XT+Z/2. )*LTAP/XT 
IF (TL.LE.C.O) TL=0.0 

IF (TL.GE.LTAP) 7L=LTAP 
AC DT=DTI+2.*(RACER*T) 

AT=3. 1A159*CT**2/A. 

CALL AREAS 

IF (ABS ( 2 Vv ) .GT.C.O) GO TC 2C 
IF(SUMAB.LE.O.C) GO TC 31 
X=(ABPCRT4ABSLC T )/SUMAB 

90 MNCZ^AT^X/APNCZ 3 * (2.* ( 1 . + B0T/2 . *MN 1 *HN 1 ) / TCP ) **ZA P 
IF (ABS (MNOZ-MNl) .LF. 0.002) GO TC 2 
MM=M!NCZ 
GC TO 90 

2 VNCZ=GAM4CSTAR*MNCZ#SCRT ( ( (2./TCP)**(T0P/BCT ) >/( l . +B0T/2 . *MNCZ*MNC 
1Z) ) 

PRAT=( 1.+80T/2.*MN0Z*MNCZ)**(-GAM/B0T) 

JRCCK=AT/APNCZ 
SUMYA=CELY* (AfiP2+ABN2+ABS2 ) 

IF (Y. EG. 0.0} SLMYA=C.C 
VC=VC+SUMYA 


-k2- 
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TABLE IV-2 (Cont’d) 


IF(Y.GT.C.O) GC TO 11 

PCNCZ= (A*RHC*CSTAR*$UPA8/AT)**S 1 . / l 1 .-N ) ) * ( 1 . + (C APGAM* JROCK ) **2/2 . 
l)**(N/( i.-N) > 

PCN-PCNQZ 

PCIS=AT*PONCZ/CSTAR 

P2=PCNCZ 

PCNCz2=PCNUZ 

PNCZ=PRAT*PCNCZ 

PA = 2.^^DIS*VNCZ/ (APHEAC+APNOZ J+PNCZ 

5 PNCZ=PRAT*PCNCZ 

PHEAD=2.*KOIS*VNCZZ ( A PHEAC+APNOZ J+PNOZ 

RHEAD^A* phe ad**n 

ZIT=MCIS*X/APNCZ 
RN 1 = R H E A 0 
PHEAD2=PHEAC 

3 RNOZ=RNl-l ( RN1-A*PNCZ**N-ALPHA*Z I7**.8/ (L**.2#EXP{ B£TA*RNl*REO/ZI T 
l ) ) ) /( 1 . + AlPhA*Z IT**.8*fiETA*RHC/Z I T/ ( L** . 2*EXP (BETA*RNi*RHO/Z I T J > ) ) 

IF IABS (RM-RNCZ ) .LE.0.CC2) GO TC 4 

RNl^RNGZ 

GO TO 3 

4 AVE1=(RH£AD4RNCZ )/2. 

IF(Y.GT.O.O) GC TO 7 

RN2=RNCZ 

RH2=RH6AD 

PCN J = PCNCZ 
CPCCY=C.O 
AVE2=AVE l 

7 RN A V£ = ( RNOZ + RN 2 } /2 . 

RHAVE= {RHEAC+RE2 1/2- 

f'G£N = RH0/2.*{ ( RNOZ+RHEAO* UBPORT+ABSLOT >+2 . * A*PCNOZ**N*ABNCZ ) 

CRCY=( AVE1-AVE2) /DELY 

R8AR-(AV£l+AVE2)/2. 

GPAX=1 .0Q2*PD IS 
GPIN=0.998*MDIS 
IF(Y.GT.G.O) GC TO 12 
GPAX=I.C01*PCI S 
GPIN=C.999«N0IS 

IF(KGEN.GE.GPrN.ANC.KGEN.LE.GMAX) GO TO 6 
PCIS^PGEN 

PCNOZ~NDIS*CSTAR/AT 
GC TO 5 

6 R6=2.*KOI$*X*L/( { APNCZ+APHEAD)*PU> 

JF(Y.LE.C.O) IxPtTElfcflOl) RE 

101 FCRPAl ( I3X, » INITIAL REYNOLDS NUMBER- •,1P£U.A) 

PCN J = PCNC Z 

17 PE^SQRT ( 2./BCT*ITOP/2.*( AE+NEl/AT )**ll./ZAP)-l.) ) 

IF{AOS(ME-PEU .LE.0.CC2 ) GC TC 9 
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TABLE IV-2 (Cont'd) 


GC TO 17 
9 CALL OUTPUT 

1C IF (Y.LE. .05*TAU ) GC TC 16 

SINKL=VC/(CAPGAMC$TAR)**2*RBAR*0PCDY/12. 

MASS=.01*PDI$ 

ANSA=Y+10.0*DELTAY 
IF (KOUNT .GT.O ) C-C TC 16 

IF(A8S(SINK1).L£. MASS. AND. ANSA. LE.ANS-XT) GC TO 18 
GC TO 16 

18 C6LY=10.*C6LTAY 
GC TO 55 

16 CELY=CELTAY 

55 CELTAT=2.*DELY/ (RHAVE+RNAVE ) 

Z=Z+DELTAT*(RNAVE-RHAVE) 

Y=Y+DELY 
T = T + OELT AT 
SUM2=SUMAB 
RN2=RNCZ 
RH2=RHE AC 
AVE2 = AVE 1 
GC TO 1 

11 MCI$^AT*PCNCZ/CSTAR 
GC TO 5 

12 DPC0Y=(PHEAC2+FCNQZ2)/( RNAVE+RHAVE )*CRDY+ ( PHEAD2 + PCNOZ2 ) / ( ( ABP2+AE 
1N2+ABS2 ) *2. )*DACY 

IF (ABS (CPCDY) .GE.DPOUT .OR.Y.GE.XCUT) GO TC 25 

SINKl=VC/(CAPGAP*CSTAR ]**2*RBAR*CPCDY/12.+{ PHEAD2+ PCNCZ2 ) /2.*I RNAV 
1E+RHAVEI/2.* t AEP2+ABN2+ABS2)/( 12.* ( CST AR*CAPGAM ) **2 ) 

STUFF = FGEN'-SINK1 

IFISTUFF.GE. GRIN. ANC. STUFF. LE.GMAX) GC TO 1A 
MC I S=STUFF 
PCNOZ=MDIS*C$TAR/AT 
GC TO 5 
1A P 1 = PCNCZ 
PCN J= PCNOZ 
PCNCZ2M PI + P2 )/2. 

P2=PCNCZ 
P3=PHE AD 

PH£AD2=( P3 + PA )/2. 

PA = PHE AD 

ANS=TAU-A8S(Z/2. ) 

IF(Y.LI.ANS) CALL OUTPUT 
IF(Y.LT.ANS) GC TO 10 

19 ZH=Z 
YV> = Y 

sumba=sumab 

P 1=PCNCZ 
RF2-KHE AD 


/ 
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RN2=RNCZ 
RAVE=AVE 1 
ABPA I N=SLP AB 
AB TC=C .0 
WRITE (6, 51) 

51 FCRPAT(A8X t **$****+*******#********• * / , A 8X » * ** ** T A 1 1 OFF BEGINS*** 
l* , t/»48Xt*’ < "************** lM * 1, ‘**** , t//) 

20 ANS2=T AU + ABS ( ZW/2. ) 

K CUNT = KOUM + I 

CEIYW=CELTAY 

CY2=0ELYW 

IF(ZW) 32*32,33 

32 IF (Y.LT.ANS2.ANC.A8SC2W) .GT.DY2) GC TC 211 
5UPA8=ABFAIN 

GC TO 31 

211 SUFCY=SUPDY+CELYW 

StiPAB*{ L » + SU^CY/ZW-CELYW/(2.*ZW) ) *ABTC- ( SUPCY/ZW-DELYW/ { 2 ,*ZW ) >*AB 
IPAIN 
GG TO 31 

33 IF(Y.LT. ANS2.ANC.ZW.6T .CY2 ) GC TC 21 
SUPA8-A8TC 

GC TO 31 

21 SUPUY-SUPCY+DELYW 

SUPAB® ( 1 ,-SUPCY/ZW+DELYW/l 2.*ZW) ) * ABM A I N+ ( SURE Y/ ZW-DEl YW/ I 2 . *Z W ) >* 
I.A0TC 

31 IF(SUNA8«LE.0.C ) PDNCZ* PONCZ/2 . 
rF { SUPAB . LE .O.C ) GC TC 25 

PCNCZ*( A*RHC*C$TAR*SuPAB/AT >**< l./l l.-N ) ) 

FCIS=AT*PCNCZ/CSTAR 
AeAVE=(SUPAB + SLP8A )/2. 

SUPYA=DELY*ABAVE 

VC=VC+SUPYA 

DADY=(SUPA6-SUPEA)/C6LY 
PEARMPl+PCNGZ ) / 2 . 

SUP8A=$UPAB 

22 CPCCY=PBAR/( l.-N)*l./ABAVE*DADY 
IF(PGNQZ.Lt.3C.C) GC TC 25 
RNCZ=A*PCNGZ**N 

rhead^kncz 

R8AR=(RHEAC+RAV£ )/2. 

PGEN=RH0*{RNG2+PHEAC )/2.*$LPA8 
GPAX=1.002*PCIS 
GP IN=0.998*P0IS 

S INK1=VC/ (CAPGAP*CSTAR )**2*RBAR*DPCCY/12 . 

STUFF-PGEN-SINRl 

IF (STUFF. GE.GPIN.ANC. STUFF. LE.GPAX) GC TO 23 
PC I S® S TUFF 

PONCZ=PONJ*CPCCY*OELY 
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IF (PONCZ.LL.C.C) PONOZ^Q .0 
PBAR- ( Pl+PCNCZ )/2. 

GC TO 22 

23 RFAVE= (RH2+RHEAC )/2. 

RNAVE= (RN2+RNCZ ) /2. 

RH2=RHEA0 
RN 2=RNCZ 
PHE AD=PGNCZ 
RA VE=RHEAO 
PI =PONCZ 
PCN J=PCNCZ 

IFIABS (DPCDY) .GE.OPGLT) GO TO 25 
IF (Y.GE.XCUT) GC TC 25 
CALL OUTPUT 
GC TO 10 
25 SLiPAB^O • 0 
Rh£ AD=C • 0 
RNCZ=RHE AD 
PHEAD=PONOZ 
WR I T£ ( 6, 318 ) 

318 FORPAT (44X, *****8EGIN HALF SECOND TRACE**** ' , // ) 

CALL CUT PUT 
TIFE=T 
CEITAT=. 5 
T I P= T I FE + 5 . 

PHT=PHEAC 

PCNT=PCNOZ 

SG=0.0 

29 T = T + DELT AT 

PHEAD=PHT/£XP tCAPGAP**2*AT*CSTAR/VCF*l T-TIP!E)*12. ) 

PCNCZ^PHEAO 

KO I S^PCNCZ * AT /C 5 TAR 

Y=Y+.5*RHEAC 

CALL OUTPUT 

20 IFIT.LT. TIM.ANC.PHEAD.GE.3C.0) GC TC 29 
IOC WP1=G*SUPPT 

HP2=RHC*(VC-VCI)*G 
ViP= { WP 1 + WP2 ) /2 • 
rsp=iTCT/HP 
ISPVAC=I TVAC/kP 

HR ITEt 6» 102 ) fcPl ,WP2, WP,PHPAX, ISP, ISP VAC, HOT, IT VAC 
102 FCRPATt 1 3X, *WP1= • , 1PE 1 1 . 4 , / , 1 3X , * WP2= • , 1PE 1 1 .4 ,/ , 13X , * WP^ , *1PE1 
11.4*/ 1 13 X f * PHP AX = • , 1PE11.4,/, 13Xt ' ISP= • , 1PE 1 1 . 4 , / , 13X , * I SPVAC= * 
2*lPEil.4*/,i3X* , ITCT- • , 1PE 1 1 . 4 , / , 1 3X , • ITVAC= S1PE11.4) 

RE AD ( 5,600) P I F K , DTEP P , S I GP AP , S I GP AS ,N 1 ,N2 , SYCNOP , DCC , PS IC »DE LC , LC 
lC,KSEG t HCN,SYNNCP,PSIS,PSIA,Kl,K2,PSI INS , CEL I NS , K6H , KEN ,DL INER , TAU 
2L.KA 

C ********************************************************************* 


-U6- 



TABLE IV-2 (Cont'd) 


C * READ IN BASIC PROPERTIES REQUIRED FOR WEIGHT CALCULATIONS * 
C * * 
C * PIPK IS THE TEMPERATURE SENSITIVITY COEFFICIENT OF PRESSURE * 
C * AT CONSTANT K * 
C * DTEMP IS THE MAX EXPECTED INCREASE IN TEMPERATURE ABOVE * 
C * CONDITIONS UNDER WHICH MAIN TRACE WAS CALCULATED IN * 
C * DEGREES FAHRENHEIT * 
C * SIGKAP IS THE VARIATION IN PHMAX * 
C ♦ SIGMAS IS THE VARIATION IN CASE MATERIAL YIELD STRENGTH * 
C * M IS THE NLMBER OF STANDARD DEVIATIONS IN PHMAX TC EE USED * 
C * AS A OASIS FOR DESIGN * 
C * N2 IS THE NUMBER OF STANDARD DEVIATIONS IN SY TO BE USED AS * 
C * A BASIS FOR DESIGN * 
C * SYCNOM IS THE NOMINAL YIELD STRENGTH OF THE CASE MATERIAL * 
C * IN LBS/ INCH * 
C * OCC IS THE ESTIMATED MEAN CIAMETER OR THE CASE IN INCHES * 
C * PSIC IS THE SAFETY FACTOR ON THE CASE THICKNESS * 


60C FORMAT (7X»F5.A»9X»F6*2f IOX » F5 . 3 » 1 OX , F5. 3 , 6X » F 5 .2 » / ♦ 5X , F5 . 2 » 1 CX ,F 1C 

K2,7X,F6.2,8X,F5.2,eX t F5.3,/,6X,F7.2,8X,F3.0,7X,F4. I , 10X „ F 10 . 2 , 8X , 

2F5.2»/,7X,F5.2t6X,F6.A»6X, F6.A, 10X,F5*2, ICX« F6 .A , / ,6X» Ft .A,7X,F6.A 
3.10X,F6.A,8X,F6.A»6X»F7.2> 

C * DELC IS THE SPECIFIC WEIGHT OF THE CASE MATERIAL IN LBS/IN**3 * 

r * LCC IS THE LENGTH OF THE CYLINDRICAL PORTION OF THE CASE * 

C * INCLUDING FCRWARD AND AFT SEGMENTS IN INCHES * 

C * NScG IS THE NUMBER OF CASE SEGMENTS * 

C * HCN IS THE AXIAL LENGTH OF THE NOZZLE CLOSURE IN INCHES * 

C * SYNNOM IS THE NOMINAL YIELD STRENGTH Of THE NOZZLE MATERIAL * 

C * IN LBS/INCH * 

€ * PSIS IS the SAFETY FACTOR ON THE NOZZLE STRUCTURAL MATERIAL * 

C * PSIA IS THE SAFETY FACTOR ON THE NOZZLE ABLATIVE MATERIAL * 

C * K 1 AND K2 ARE EMPIRICAL CONSTANTS IN THE NOZZLE WT. EQUATION * 

C * P S 1 INS IS THE SAFETY FACTOR ON NOZZLE INSULATION * 

C * CE L INS IS THE SPECIFIC WEIGHT C-F THE INSULATION IN LB$/1N**3 * 

C * KEF IS THE EROSION RATE OF INSULATION TAKEN CONSTANT * 

C * EVERYWHERE EXCEPT AT THE NOZZLE CLOSURE IN IN/SEC * 

C * KEN IS THE ERCSICN RATE OF INSULATION AT THE NOZZLE CLOSURE * 

C * IN IN/SEC * 

C * CLINER IS THE SPECIFIC WEIGHT Cf THE LINER IN LBS/IN**3 

C * TAUL IS THE THICKNESS OF THE LINER IN INCHES 

C * WA IS ANY ADDITIONAL WEIGHT NCT CONSIDERED ELSEWHERE IN LBS * 

r ifr*****************************^************************************* 

WRITE (6, 610) PIPK,DTEMP,SIGMAP,SIGMAS,N1,N2,SYCNCM,CCC,PSIC,DELC,L 
1 CCtNSEG,HCN,SYNNCM, PS IS, PS I A t Kl ,K2,PS II NS, DEL INS »KEH, KEN, CLINER, TA 
2UL » WA 

6 1 C FORMAT {//,2CX, • INERT WEIGHT INPUTS' »/* 13X r 'PIPK- « , 1PE 1 1 . A ♦ / , 1 3 X » • 
1CTEMP= » ,lPfcll.A T /,13X, *SIGMAP= * ♦ 1 P E 1 1 . 4 * / , 1 3X , ' S I GMA S= »,IPE11.A 
2 » / 1 1 3X » * Nl= • , IPEi 1.4,/, 13X, »N2= • , 1PE U . A , / , 13X , « SYCNCM= SlPfclL. 
3A » / , 1 3X , ' CCC= , ,1PE11.A,/,13X,*PSIC= * , 1PE 11 .4 ,/ » 13X , • CELC* SIPE1 


* * 
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A 1 • A » / » 13X , • LCC= ' , 1PE1 l.A, /, L3X, *NSEG= * , IPE l l .A , / , 13X , *HCN= *,IPE 
51l.A,/,l3X, 'SYNNCM* ' » IPE 1 1 . A * / , 1 3X , • PS I S* ' « 1 PE 11 . A , / , 13X , • PS I A = 
6'»IPE11.A,/ »13X»'K1- • , 1PE11.A,/,13X, 'K2^ « , IPE 1 1 . A * / # 1 3X , • PS II NS= 

7 *,lPEll.A,/,i3X,«CELINS= ■ , 1PE11 .A , / , 1 3X , • KEF= • , 1PE1 1 .A, /, 13X » • K 

8 EN= • , 1PE11.A,/ , 13X, »DLINER= ♦ , IPE 11 . A , / , 1 3X , • TAUL= • , IPE 1 1 .A , / , 13 
9X,*WA= 1 , IPE 1 1 • A ) 

PMECP=PHMAX*( l.*Nl*SIGMAP)*EXP(PIPK*DT£MP) 

SYC=SYCNOM* ( 1 .-N2*SIGMAS ) 

TAUCC=PS IC*PMECP*DCC/( 2.*$YC) 

WCC=3 • 1A159#TACCC*CCC*CELC*LCC* ( 1 . + ( NSEG- 1 . ) * I AO .*TAUCC /LCC ) ) 
TAUCD=lAUCC/2. 

WCH=2.5*3. 1A155/2.*CCC**2*TAUCC*DELC 
WCN=A.5*3. 1A155/2.*DCC*HCN*TAUCD*DELC 
WC=WCC+WCH+ViCN 
EPS IL=AE/A T 

DT=2.*SCRT ( AT/3* IA159) 

WN=Kl*DT**2/( 1 • + .5*S IM ALFAN ) )*( ( EPS I L-SGRT ( EPSIL ) )*PMEOP*DT*PSIS/ 
1SYNN0M+K2*T*PSIA ) 

WINS = T*PSI lNS*CELINS*CCC*3.1A159*tKEH*(DCC*.AC+( S+NS )*TAU/2.+0. 15/ 
1PSIINS*(LCC-TAI*(S+NS)) ) +KEN* . 80*FCN ) 
KL=TAUL*OLINER*3.1A155*DCC*(OCC/2.+LCC+HCN) 

W'M=WC + WN + W I N S + W L + WA + WP 
ZETAM=WP/WM 
RA T 1 0= I TCT/WM 
WRITEI6, 605) 

605 FORMAT!///, A2X,'M0TCR WEIGHT CAICULAT IONS • ) 

WRI TE 16»601 ) PMEOP, TAUCC,WC,WN,WlNS,WL,WM,ZETAM,RATIO 
601 FORMAT { 13X, 'MAX EXPECTED PRESSUR£= • , IPE 11 .A, / , 13X , * CYL INCR ICAL CA 
1SE THICKNESS^ • , 1PE1 1 .A, /, 13X, • CASE WT= • , IPE 1 1. A, /, 13X , • NCZZLE WT 
2= *, 1PE11. A, /, 13X, • INSULATION WT= ', 1P611. A,/* 13X, 'LINER WT= SIPE 
311.A,/,13X, • TCTAL POTCR WT= • , IPE U.A ,/ , 13X , • ZET AM= * , IPE 1 1 . A, / , 1 3 
AX, 'RATIO CF I TCT TO WM = 1PE11. A) 

901 CONTINUE 
STOP 
END 
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C 

C 

c 

c 

c 


c 

c 

c 


SUBROUTINE AREAS 

* 44 #***«*«* 4 ** 4 ** 4 * 4 ********* 4 ******* 4 ********** 4 ************** 4 ***** 

* SUBROUTINE AREAS CALCUUAT6S BURNING AREAS AND PORT AREAS FOR * 

* CIRCULAR PERFORATED (C.P.) GRAINS AND STAR GRAINS OR FOR A * 

* COMBINATION CF C.P. AND STAR GRAINS * 

INTEGER STAR t GRA IN, ORDER 

REAL MGEN,MCIS,MN0Z,MM » JRCCK , N » L , ME 1 ,ME, ISP, I TO T , MU ,MASS » ISPVAC 

REAL LGCI ,LGNI ,NS,NN,NP,LGSI ,NT , L T P , LGC * LS » LF 

REAL M2, MDBAR, ISP2, ITVAC 

C0MM0N/CCNSTI/2R ,AE, AT , THETA, ALFAN,G 

CGMMCN/CCNST3/S ,NS 

COM MON /VAR I Al/ V , T , CE L Y , CEL T AT , PCNC Z , PHE AD , RNCZ ,RHE AC » SOM A8 , PHM AX 

C0MMQN/VARIA2/ABP0RT , A BS LO T » ABNCZ » APHEAC,APNCZ ,D ADY, AEP2 , A0N2 , ABS2 

COM MON /VAR I A3 / I TOT, I T VAC , JRCCK , I S P , I SPV AC ♦ MD I S , MNC Z , SG , SUMMT 

CCMM0N/VARIA4/RNT,RHT,SUM2,R1,R2,R3,RHAVE,RNAV£,RBAR,YB,K0LNT,TL 

COM MON/ VAR IA5/AEMAIN»ABT0,SUMDY 

A8PC-C .0 

ABNC=0 . 0 

AESC=C .0 

AB PS=C .0 

ABNS=0 .0 

ABS S=0 .0 

DA8T=C.O 

BBP^O.O 

B BN=0 . 0 

BBS^O.C 

SG=Q.C 

ANUM=3. 14153/4 . 

RN'T = RNT + RNOZ*DELTAT 
RFT=RHT+RHtAC*CELTAT 

1 IF (ABS ( Z to ) .LE.C.C) K=C 
IF (ABS (Zto).GT.C.O) K=I 
YB = Y 

IF(K.EC.l) Y=YE-SUMCY/2. 

2 IFIK.EC.2) Y=YB+aBSI ZKJ/2.-SUM0Y/2. 

IF(Y.LE.C.O) REAC(5,5CC) I NPUT , GR A I N , STAR , NT , CROER 


* READ THE TYPE OF INPUT FOR THE PROGRAM AND THE BASIC GRAIN * 

* CCNFIGURAT ICN AND ARRANGEMENT * 

* VALUES FCR INPUT ARE * 

* 1 FOR CNLY TABULAR INPUT * 

* 2 FOR CNLY EOLATION INPUTS (EQUATIONS ARE BUILT * 

* INTO THE SUBROUTINES * 

* 3 FCR A COMBINATION CP I ANC 2 * 

* VALLES FCR GRAIN ARE * 

* L FCR STRAIGHT C.P. GRAIN * 

* 2 FCR STRAIGHT STAR GRAIN * 
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C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


* 3 FOR COMBINATION OF C.P. ANC STAR GRAINS * 

* VALLES FOR STAR ARE * 

* 0 FOR STRAIGHT C.P. GRAIN * 

* I FOR STANDARD STAR * 

* 2 FOR TRLNCATEO STAR * 

* VALUES FCR NT ARE * 

* C IF THERE ARE NO TERMINATION PORTS * 

* X WHERE X IS THE NUMBER CF TERMINATION PORTS * 

* VALUES CF ORDER ESTABLISH HOW A COMBINATION C.P. AND STAR * 

* GRAIN IS ARRANGED * 


* 1 IF DESIGN IS STAR AT HEAD END AND C.P. AT NOZZLE * 

* 2 IF DESIGN IS C.P. AT H6AC ENC AND C.P. AT NOZZLE * 

* 3 IF DESIGN IS C.P. AT HEAD END AND STAR AT NCZZLE * 

* A IF DESIGN IS STAR AT HEAD END AND STAR AT NOZZLE * 

* ***NCTE*** IF GRA l N = 1 1 VALUE OF CRCER MUST BE 2 * 

* ***NCTE*** IF GRAIN=2, VALUE OF ORDER MUST BE A * 

500 FORMAT (9X,I2»9X,I2*8X t l2f6X*F4.0,9X, 12) 

IF (Y.LE.C.O) WRITE(6,607) 

607 FORMAT {//»20X» ’GRAIN CONFIGURATION*) 

IFIY.LE.O.O) WRITE(6,6C0> I NPUT ,GRAI N , STAR ,NT , OROER 
600 FORMAT ( L3X» * INPLT= * , 1 2 , / , 13X , • GRA IN= • f I 2, /, 13X , * STAR* %I2,/,13X 
l » * NT= * * F4 . C » / tl3X»*CRCER= « t I2,//) 

I F ( INPUT .EQ.2 ) GO TO 12 
IFIY.LE.O.O) GO TO 6 
IFtYT.LE .Y.ANC.K.LT.2) GO TO 8 
g DENGM* YT-YT 2 

SLOPE 1= ( AB PK-AEFK2 ) /DENOM 

SLCPE2=( ABSK-AESK2)/DENCM 

SLOPE 3=(ABNK-AeNK2)/D£NCM 

SLCP£4*( APHK-AFFK2) /DENGM 

SLOPES*! APNK-AFNK2)/DENCM 

B1=A8PK-SLUPE1*YT 

82=ABSK-SL0PE2*YT 

B3=ABNK-SLCPE3*YT 

B4=APHK-SLCPE4*YT 

B5=APNK-SLOPE5*YT 

ABPT=SLOPEl«Y+ei 

ABST=SL0PE2*Y+B2 

A8NT*SL0PE3*Y+E3 

APHT=SL0PE4*Y+e4 

APNT=SL0P£5*Y+E5 

IFIINPLT.EQ.3) GO TO 3 

GC TO 52 

6 READI5»507) YT , A BPK f AB SK , ABNK , APhK , APNK 

444*444*444‘4*#44444**4444**4*44«4***4***4**+44****$*4*******4**4*’»'l‘*4 

* READ IN TABULAR VALLES FOR Y=0.0 (NOT REQUIRED IF INPUT=2) * 

* * 
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c 

* 

ABPK 

IS 

THE 

BURNING AREA 

IN 

THE PORT IN IN**2 

* 

c 

* 

ABSK 

IS 

THE 

BURNING AREA 

IN 

THE SLOTS IN IN**2 

♦ 

c 

* 

ABNK 

IS 

THE 

BURNING AREA 

IN 

THE NOZZLE END IN IN**2 

* 

c 

♦ 

APHK 

IS 

THE 

PORT AREA AT 

THE 

HEAD END IN IN**2 

* 

c 

* 

APNK 

IS 

THE 

PORT AREA AT 

THE 

NCZZLE ENC IN IN**2 

* 


C *«« 44 44**** « 4444***4* 4**4 ****444** ****** 44444444(4444444444444 *4444*44 
507 FORMAT l6X,Ffc.2»10X,eil.*,iCX,E11.4,8X,£11.4,/»22X,Ell.<t t 9X,E11.4) 
KR I TE ( 6 » 6 10 ) 

610 FORMAT (/ »13X, • TABULAR VALUES FOR VT ECUAL ZERO READ IN') 

A8PT=A8PK 

ABST=A8SK 
ABNT=A6NK 
APHT=APHK 
APNT=APNK 
YT 2 = YT 

IF (INPUT. EG.3 ) GC TC 3 
GC TO 52 
8 YT2-YT 

ARPK2=ABPK 
ABNK2=ABNK 
APSK2= AQSK 
APHK2-APHK 
APNK2=APNK 

RE AD I 5, 5 05 ) YT , A BPK , ABS K , A SNK , APHK , APNK 
C 4 **** 444 * 44444 * 4444 * 4 * 4 ** 44 * 444444 * 44 *** 4444 * 4 * ******* 44 * 4 ***** 444*44 

C * READ IN TABULAR VALUES FOR Y=Y ( NCT RECU I R ED FOR INPUT=2) * 

c * 4 * 44444 * 444 * 4 *** 4 * 4 **** 4 ******** 4 *** 44 ** 4 ******* 4 ***#****** 4 *** 4*444 

505 FCRMAT{6X,F6.2,10X»ElI.4,iCX,EIK4,8X,Ell.4*/,22X f en.4,9X,EU.4) 
WRITE<6,6in YT 

611 FCRMAT(/,13X, 'TABULAR VALUES FOR YT= ',F7.3 9 ' READ IN') 

GO TO 9 

12 AfiPT=C«0 
A0NT=C.O 
ABST=0.0 

3 IF (GRAIN. Nt. 2) GC TG A 
AEPC=C.O 

A ENC = 0 .0 
AB$C=0.0 
GO TO 7 

4 IF(Y.LE.O.O) READ(5,5C1) DC ,DI , DELD1 * S, THETAG »LGC I » LGN! , THE TCN, THE 
ITCH 

*<****4**4****44444**4*****44444****4*******4**44*44**4******44****** 


* READ IN eASIC GEOMETRY FOR C.P. GRAIN (NCT RECUIREC FOR * 

* STRAIGHT STAR GRAIN) * 

* DO IS THE AVERAGE CUTSICE INITIAL GRAIN DIAMETER IN INCHES * 

* Cl IS THE AVERAGE INITIAL INTERNAL GRAIN DIAMETER IN INCHES * 

* CELCI IS THE DIFFERENCE BETWEEN THE INITIAL INTERNAL GRAIN * 

* OIAMETER AT THE NOZZLE ENC AND Cl IN INCHES * 


• 51 - 



GE 13 


TABLE IV-2 (Cont'd) 


C * SIS THE NUHtER CF FLAT BURNING SLOT SIDES (NOT INCLUDING ♦ 
C * THE NOZZLE ENC) * 
C * THETAG IS THE ANGLE THE NOZZLE END OF THE GRAIN HAKES WITH ♦ 
C * THE NCTCR AXIS IN RADIANS * 
C * LGCI IS THE INITIAL IOTAL LENGTH OF THE CIRCULAR PERFORATION * 
C * IN INCHES * 
C * LGN I IS THE INITIAL SLANT LENGTH OF THE BURNING CONICAL * 
C * GRAIN AT THE NCZZLE END IN INCHES * 
C * THE TC N IS THE CONTRACT ICN ANGLE CF THE BCNCED GRAIN IN RAC. * 
C * ThETCH IS THE CONTRACTION ANGLE AT THE HEAC ENC IN RADIANS * 


501 FORPAT <5X T F7.3,6X,F7.3,9X,F7.3,5X,F4.C»9X,F7.5,/t7X,F7.2,7X,F6.2,9 
1X,F7.5,9X,F7.5) 

IF(Y.LE.C.O) WRITE (6 ,601) CO, 01 ,CELOl , S , THETAG ,LGC I ,LGNI , THETCN , TH 
1ETCH 

601 FOR PAT (20X, *C . F • GRAIN GE0RETRY*,/,13X, *DC= * , F7 .3 , /» 13X, * D 1= ' ,F7 
1.3,/,13X,'CELCI= ' »F7.3»/»13X,»S= *,FA. 0, /, I 3X, * THETAG = *,F7.5,/,1 
23X, • LGC I = • ,F7.2,/,13X,*LGM* «, F6 . 2 ,/, 13X ,» THETCN = * ,F7 . 5 , / , 13X f • 
3 THE! CH= * ,F7.5,//) 

BNUP=ANUP*CC**2 
T LL- TL 

I F ( ORCER . GE .3 ) TLL=C.O 
ABSC=S*ANUP*ICC**2-(CI+2.#Y)**2) 

IF ( ABSC.LE.C.O ) ABSC=C.C 
IF (2.+Y+CI.GT.CC ) GO TC 100 
IF ( THETAG. Lfc.C .C8727 ) GO TC 101 

ABPC=3.i4159*(CI+2.*Y)*{ LGC I-(S+TANtTHETAG/2.) ) *Y-Y*COTAN I THETCH )~ 
1TLL) 

I F ( AB PC . LE . C • 0 ) AB PC-C . 0 

ABNC=3. 14i59*< LG N I - Y *CC T AN ( THET AG+ THETCN > - Y*T AN ( THETAG/2. I }* IDI+DE 
1LCI+Y4LGNI*SIN(THETAG)+Y/SIN ( THETAG+THETCN)*S INI THETCN) ) 

IF ( ABNC.LE.C.C ) ABNC=C.O 
GC TC 5 
IOC ABNC=C .0 
ABPC= C . 0 
GC TO 5 

101 LGC=LGCI~Y*(CGTAN( THETCN )+COTAN( THETCH) ) 
A8PC=3*14159*(CI+2.*Y)*(LGC“S*Y-TLL ) 

IF ( ABPC.LE .C.C ) AB PC=C • C 
ABNC=C.O 

5 APHT=ANUM*ICI42.*RHT)**2 
I F { APHT .GE. 8NLP ) APHT=BNUH 
IFIK.LT.2) APHT 1=APHT 
APNT=ANUP*(Cl+CELDI+2.*RNT>**2 
IF ( APN r . GE . BNLP ) APN T = BNUH 
I F ( GRA IN *NE • 1 ) GC TC 7 
ABPS=C.O 
AESS=C.O 


y 


- 52 - 


u U V l*> O O u 


ice 14 


TABLE IV-2 (Cont'd) 


ABNS=C .0 
GO TO 50 

7 IF(Y.LE.O.O) R£AD(5 f 5C2) NS » LGS ! t NP , RC , F ILL . NN 

C *4***#**********jM*********** *******+*******************#************ 


C * READ IN BASIC GEOMETRY FOR STAR GRAIN (NOT REQUIRED FOR * 
C * STRAIGHT C.P. GRAIN) * 
C * NS IS THE NUMBER OF FLAT BURNING SLOT SIDES (NOT INCLUDING * 
C * THE NOZZLE END) * 
C * LGS I IS THE INITIAL TOTAL LENGTH OF THE STAR SHAPED * 
C ♦ PERFORATED GRAIN IN INCHES * 
C * NP IS THE NUMBER OF STAR POINTS * 
C * RC IS THE AVERAGE STAR GRAIN OUTSIDE RADIUS IN INCHES * 
C * FILL IS THE FILLET RADIUS IN INCHES * 
C * NN IS THE NUMBER OF STAR NCZZLE END BURNING SURFACES * 


502 FORMAT(5X,F4.0t7X # F7.2,5X,F4.0,5X f F7.3 T SX f F7.3,5X,F4,C ) 
IF(Y.LE.O.O) WRITE <6»602) NS»LGSI*NP»RC»FlLLtNN 
602 FORMAT! 15X, 'BASIC STAR GEOMETRY 13X, »NS= • ,F4 .0 ,/ » 13X f ' LGS I- *, 
1F7.2,/,13X, 'NP= ',F4.C,/»13X, *RC* * ,F7 . 3 , / , 13X , * F I LL= *»F7.3,/,13X 
‘2,'NN* »,F4.C,//) 

FY = F I L L+ Y 

I F ( STAR . £Q . I ) GC TO 20 
IF(Y.LE.O.O) RE A D l 5 1 503 ) RP»TAUS 


* READ IN GEOMETRY PCR TRUNCATED STAR (NOT REQUIRED FOR * 

* STANDARD STAR) * 

* RP IS THE INITIAL RAOIUS OF THE TRUNCATION IN INCHES * 

* TAUS IS THE THICKNESS OF THE PROPELLANT WEB AT THE BOTTOM * 

* CF THE SLOTS IN INCHES * 


** ********* ************** ******************************************** 
503 FORMAT (5X»F7.3»7X»F7.3) 

IF(Y.LE.C.O) V»RIT£(6f603) RP,TAUS 

603 FORMAT (20X, 'TRUNCATED STAR GEOMETRY' ?/» I3Xt*RP= ’ »F7*3»/»13X » * TAUS 
1= SF7.3,//) 

THETAS=3«14159/NP 

RPY=RP+Y 

LS=RC-TAUS-FILL-RP 

RPL=RP+LS 

THETS1=THETAS-ARSIN(FY/RPY ) 

IF (THETS I.LE.C .0 ) GO TC 110 
IF(Y.LE.TAUS) GC TO 103 

THETAC=ARSIN( (RC*#2-RPL**2-FY**2)/(2.*FY*RPL) ) 

IF(THETAC.GE.C.C) GC TO 104 
IF(Y.LT.RC-RP) GO TC 105 
SG=0 , 0 
GC TO 14 

103 SG=2.*NP*(RPY*TFETSl+LS-(RPY*C0S(THETAS-THETSl)-RP 1+3. 14159/2.*FY) 
GC TO 14 


/' 
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104 SG=2.*NP*(RPY*TF£TS1+LS-(RPY*CQS(7HETAS-THETS1)-RP)+FY*THETAC) 

GC TO 14 

105 SG=2.*NP*(RPY*THETS1+SCRT(RC**2-FY**2)“SCRT(RPY**2-FY**2 ) ) 

1 A AGS=3.14 159*(RC**2-RP**2 J-NPM3. 14 1 59*F ILL**2/2 . +2 . *LS*F I LL ) 

GC TO 31 

110 THE TAF= T HE T AS 

THETAP=2.*TFETAS 
TAUWS=TAUS 
GC TO 111 

20 IF(Y.LE.O.O) REAC{ 5 » 5C4) THETAF , THETAP, TAUWS 
C 4*4 *****4*****44*44****4*4********4*******4* ************************* 


C * READ IK GECRETRY FCR STANDARC STAR (NCT REQUIRED FCR * 
C * TRUNCATED STAR) * 
C * THETAF IS THE ANGLE LOCATION OF THE FILLET CENTER IN RADIANS * 
C * THETAP IS THE ANGLE CF THE STAR POINT IN RADIANS * 
C * TAUWS IS THE WEB THICKNESS OH THE GRAIN IN INCHES * 


C 4***44*44**4**4*444444444444444*4444444444444**444*4**4*44444**4*4444 
50A FORRAT (9X,F7.5,9X,F7.5,8X,F6.3) 

IFtY.LE.O.O) W R I T£ ( 6 * 604 ) THETAF f THE TAP* TAUWS 
60A FORRAT (20X, 'STANDARD STAR GEORETRY 13X THETA F= *,F7.5,/,13X,*T 
1HE TAP- ' ,F7.5,/,13X,»TAUW$ = * »F6 • 3,// ) 

THET6S=3.l4l59/NP 
THETS 1=1.00 
ill LF=RC-TAUkS-FILL 
CNUR=(Y+FILL)/LF 
DNUR = S IN (THETAF )/S IN (THETAP/2. ) 

ENUR=(RC**2-LF**2-FY**2 )/{2.*LF*FY) 

FNUR=S lN(THETAF)/C0S(THETAP/2.) 

IF (CNLR.LE.FNUR ) GC TC 106 
I F ( Y . LE- TAUWS)GC TC 107 

SG=2. *NP*FY*(THETAF+ARSIN( SIN (THETAF > /CNUK J-ARCOS ( ENUR ) ) 

GC TO 23 

106 IF ( Y. LE. TAUWS) SG=2 .*NP* LF* ( DNUR+CNUM* ( 3 . 14 159/2 .+ THETAS- THE TAP/2 . 
1-CCTAN ( THE TAP/ 2. ) ) + THET AS-ThETAF ) 

IF(Y.LE. TAUWS) GC TC 23 

SG=2. *NP* ( F Y* A RS IN (ENLR-( THETAS- THETAP/2. ) )+LF*DNUR-FY*COT AN (THETA 
IP/2. ) ) 

GC TO 23 

107 SG=2.4NP*LF*(CNLR* ( THE TAS+AR S I N { SI N { THE T AF )/CNUM) )+THETAS-THETAF ) 
23 IF (THETS1.LE.C.C ) GC TO 14 

AGS=3. 14159*RC**2-NP*LF**2* (SIN( THETAF)* (CCS ( THETAF)- SIN ( THETAF )*C 
1CTAN( THETAP/2. ) ) +THET AS-THETAF+2 . *F I LL/LF * { S IM THETAF ) /S I N ( THE T AP/ 
22. >+THETAS-THETAF+FILL/{2.*LF ) * ( 3 . 14 159/2 .+THETA S-THETAP/2 .-COT AN ( 
3THETAP/2. ) ) ) ) 

31 IF (SG.LE.C.O) SG=0 • C 
IF(K.EC.O) SGN=SG 
IF(K.LE.l) $GF*SG 
IF(K.£G.2) SGN=SG 
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IF(Y.LE.C.O) SG2-SG 

IF(K.EC.2> GO fC 37 

RAVEDT=Rl+(SG+SG2)/2.*R8AR*DELTAT 

RNDT=R2+(SG+SG2)/2.*RNAVE*DELTAT 

RHDT = R3+($G-*SG2)/2.*RHAV£*CELTAT 

Rl=RAVEDr 

R 2 = R N C T 

R3=RHCT 

SG2-SG 

GC TO 38 

37 I F { KOUNT • NE . 1 ) GC TC 39 
SG 3= SG 

RA = R1 
R5 = R2 
R6 = R3 

39 RAVEDT=RA+( $G+$G3)/2.*REAR*CEL7AT 
RN'DT = R5 + l SG+SG3 )/2.*RNAVE*CELTAT 
RHCT=R6+ ISG+SG3 )/2.*RHAVE*DELTAT 
RA=RAVEDT 
R5=RNCT 
R6 = RHC T 
SG 3-SG 

38 A8SS= ( AGS-RAVECT 1*NS 
IF(ABSS.LE.C.O) ABSS=C.C 
IF(SG.Lb.C.O) ABSS=C.C 
A8NS=(AGS-RNDT ) *NN 

IF (ABNS.LE.C.C) ABNS=0.0 
IFISG.Lt.C.C) ABNS=G.C 

IF (ORDER. LE. 2) ABP$= { LGS I- Y* ( NS+NN ) )*SG 
IF(CRCER.Lfc.2) GC TC 36 
ABPS= (LGSI-TL-Y* (NS+NN) )*SG 
36 APHS=3.1A159*RC**2-AGS+RHDT 

TF ( APHS.GE.3. 14159*RC**2) APHS=3 . 14159*RC**2 
IF (SG.LE.O.O) APHS=3.I4159*RC**2 
APNS=3.I 4159*RC**2-AG$+RND7 
IF(K.LT.2> APHS1=APHS 

If {APNS.GE.3. 14159*RC**2) A PN S=3 . 14 1 59*RC**2 
5G IF (NT.EG.C.O) GC TC 52 

If(Y.lfc.O.O) R£ A C ( 5 , 5C6 ) L TP, CTP, THETTP , TAUEFF 

C ***$*****#********J»#**$**$*******iC!*#**j<>******#***$**##*******jC!*jS*j»##* 


C * READ IN G6CM E TRY ASSCCIATEO MTh TERMINATION PORTS (NOT * 

* REQLIREC IF NT=0) * 

* LTP IS THE INITIAL LENGTH CF THE TERMINATION PASSAGES ♦ 

* IN INCHES * 

* DTP IS THE INITIAL DIAMETER CF THE TERMINATION PASSAGE * 

* IN INCHES * 

* THETTP IS THE ACUTE ANGLE BETWEEN THE AXIS OF THE PASSAGE * 

* AND THE MOTOR AXIS IN RAC I AN S ♦ 



GE 


17 


TABLE IV- 2 (Cont'd) 


C * TAOEFF IS THE ESTIMATED EFFECTIVE WEB THICKNESS AT THE * 

C * T ERR l N A T ION PCRT IN INCHES ♦ 

C *444*****«***44****4****t*>M*****4*4********************************* 

506 FORRAT (7X,F6.2,7X,F5.2, 10X.F7.5, 10X,F6.3> 

IF(Y.LE.C.O) WBITE16.6C6) LTP ,DTP, THETTP, TAU6FF 
606 F0RRA7 (20X, 'TEBR1NATICN PCRT GfcCME TRY » , / , 1 3X , * LT P= • , F6 .2 , / , l 3X, • C 
1TP = »,F5.2,/,13X,'THETTP= 1 ,F7.5,/»13X» * TAUEFF = *,F6.3,//) 

CAB T = NT*3. 14159* ( ( D TP+2 . *Y ) * ( LTP-Y/S I N ( THET TP ) )~ ( DTP+2 . *Y ) **2/4 .+ 

1 (Y+DTP/2. )*(DTF/2. )*< l.-l./S I N( THETTP) ) > 

IFtY.GE. TAUEFF ) DABT=G.O 

52 bep=o.o 

E B $ = 0 . G 
EBN=0 . C 

AEPORT=ABPT+ABPC+ABPS+CABT+EBP 
ABSLOT = A0ST + AE$C*AB$S-*EBS 
A8N0Z=ABNT+ABNC+ABNS+BBN 
SUPAB^ABPORT+AESLOT+ABNCZ 
IF(K.EC.O) GO TC 99 

IFU.ECi. 1) ABPAIN = ABPCRT + ABSLCT+ABNCZ 
K = K+ 1 

IFIK.GT.2) GO TC 69 
GC TC 2 

69 ABTC=ABPGRT+ABSLOT+ABNCZ 
99 CONTINUE 

IF(Y.GT.C.O) GC TO 70 
ABP 1 = ABPCR T 
ABN1=ABN0Z 
A8S1=AB$LCT 

70 ABP2-(ABPl+ABPCRT)/2* 

ABN2=< ABN1 + ABNCZ )/2. 

A6S2={ABSl+A8SLCT)/2. 

I F C INPUT. EG. i) GC TC 76 
GO TO (71,72,73, 74) .ORDER 

71 APHEAD=APHSi 
APNOZ^APNT 
SG=SGH 

GC TO 75 

72 APH6AC=APHU 
APNCZ^APNT 
SG=0 . 0 

IF(GRAIN.EQ.3) SG=( SGH+SGN )/2. 

GC TO 75 

73 ApHEAC=APHTl 
APNCZ= APNS 
SG=SGN 

GC TO 75 

74 APHEAD=APFS1 
APNCZ= APNS 
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TABLE IV-2 (Cont'd) 


SG = SGN 
GC TO 75 

76 APHEAC=APHT 
APNCZ=APNT 

75 Y = YB 

D!FF=SLNAB-SUN2 
DACY = C IFF/DELY 
ABP 1=A8P0R T 
A8N1=ABNCZ 
AB S 1 = A8S 10 [ 

IF(ZW.GE.Q.C) GC TC 77 
ABP1=ABMAIN 
ABFA lN=ABTO 
AB TC = ABP 1 

77 RETURN 
END 
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TABLE IV-2 (Cont’d) 

SUBROUTINE CUT Ft T 

C ^* 4 ******************** tf ****** ********* ****************************** 


C * SUBROUTINE output calculates basic performance parameters * 
C * AND prints THEM CUT * 
C * (WEIGHT CALCULATIONS ARE PERFORMED IN THE MAIN PROGRAM) * 
C * T IS THE TIME IN SECS * 
C * Y IS THE DISTANCE BURNED IN INCHES ♦ 
C * RNCZ IS THE NOZZLE END BURNING RATE IN INCHES/SEC * 
C * RHEAO IS THE HEAD ENO BURNING RATE IN INCHES/SEC * 
C * PCNOZ IS THE STAGNATION PRESSURE AT THE NOZZLE END IN PSIA ♦ 
C * PHEAD IS THE PRESSURE AT THE HE AC END CF THE GRAIN IN PSIA * 
C * PTAR IS THE PORT TC THROAT AREA RATIO * 
C * MNCZ IS THE MACH NUMBER AT THE NOZZLE END CF THE GRAIN * 
C * PATH IS THE ATMOSPHERIC PRESSURE AT ALTITUCE IN PSIA * 
C * SUMAB IS THE TOTAL BURNING AREA CF PROPELLANT IN [N**2 * 
C * SG IS THE BURNING PERIMETER IN INCHES CF THE STAR SEGMENT * 
C * (IF ANY) * 
C * CFVAC IS THE VACUUM THRUST COEFFICIENT * 
C * FVAC IS THE VACUUM THRUST IN LBS * 
C * F IS THE THRUST IN LBS * 
C * F IS THE THRUST IN LBS * 


C $#$*:(i!Mi*#*«j»>**:Mi<'******$$**$#****$***« ****#*♦#******************* + ♦** 
REAL MGEN»MClS,MNOZ»MM,JRCCK,N,L»MEl,ME,ISP,ITOT,MU*MASS»ISPVAC 
REAL M2, MDBAR, ISP2, ITVAC 

com mon /const i/zw , ao, at,theta,alfan,g 

C0MMON/CCNSr2/CAPGAM,ME,8OT,ZErAF,TB,HB,GAM 

COMMDN/VARI A 1/Y,T, CEL Y , CEL TAT, PCNCZ, PHEAD, RNCZ, RHEAD, SUMAB, PHMAX 
COM MON /VAR I A3/ I TOT, ITVAC ,JRCCK, I SP , I SPVAC ,MDI S ,MNOZ, SG.SUMMT 
YB=*Y 

IF(Y.LE.O.O) M2=MDIS 
M.DBAR= (M2 + MCIS )/2. 

SUMMT-SUMMT+MCBAR*D£LTAT 
PTAR = i ./ JRCCK 

PR£S=(I.+BOr/2.*ME*ME)**<-GAM/80T) 

ALT = HB*(T/TB}**(7*/3. ) 

PATM^IA. 696/EXP ( C . 43 1 C3E-04*AL T ) 

CF=CAPGAM*SCRT(2.*GAM/BCT*< U-PRES**(80T/GAM) ) ) + AE/AT* ( PRES-PAIM/P 
1CNCZ ) 

CFVAC=CF+AE/AT*PATM/PCNCZ 

F=ZETAF*CCS(THETA)*PCNCZ*AT*( ( 1 . +COS ( ALFAN ) ) /2 .*CF+ ( 1 .-COS ( ALFAN ) ) 
I/2.*A£/AT*IPR£S-FATM/PCN0Z) ) 

IF(F.LE.C.O) F=C.O 
IF(Y.LE.C.O) f 2 = F 
FBAR=(F+F2)/2. 

FVAC=ZETAF*CCS( THETA )*PCNGZ*AT*( ( I . +CCS ( A LF AN ) ) / 2 . *CFV AC* ( I . -CCS ( A 
ILF AN) )/2.*AE/AT*PRES) 

IF(Y.LE.O.O) FV2=F VAC 
FVBAR=(FV2+FVAC) /2. 
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TABLE IV-2 (Coat'd) 


I$P=F/(MCIS*G) 

ISPVAC=FVAC/ (PC IS*G ) 

IT0T=1 TOT+FeAR*CELTAT 
I T VAC= ITVAC+FVEAR*DELTAT 
P2 = PD I S 
F2 = F 

FV2=FVAC 

IF(PHEAD.GT.PHPAX) PHR AX=PbEAC 

HR I TE ( 6 » 1 > T f YE|RNOZ»RFIEAD, PCNOZ f PHEAC » PT AR , RNCZ » SUPAB » SG » PA TR 1 ,CF V 
1 AC , FVAC » F 

1 FCRMATU3X, *TIPE= • ,F5. 1 , 12X » ' Y= • ,F6.3, / , 13X, *RNCZ= • » 1PE 1 1 .4 . ' 
1RFEAD= 'tlPEll.4,* PCNOZ= ',1PE11.4,» PHEAC= • , 1PE 1 1 . 4 , / 1 1 3X , • P T 
2AR= ' t 1PE 11 .4 , 1 RNCZ = MPEU.4,* SUKA6= •,1PE11.4,> SG= •. 

31PE11.4, /,13X, 'PATP* 1 » 1PE 1 1 * 4 » • CFVAC= SlPEli.4,' FVAC= *,1PE 
4H.4,' F= **1PE11.4,//) 

RETURN 

END 
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Figure IV-1 (Cont’d) 
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Figure IV-1 (Cont'd) 
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Figure IV-1 (Cont’d) 
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Figure IV-1 (Cont'd) 
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Figure IV-1 (Cont’d) 
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Figure IV-2 (Cont'd) 
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Figure IV-2 (Cont'd) 











-71“ 



= THETAS 

s 

L f = LF 

9 THETS1 
sx 

L = LS 
s 


S G " SG 


6 = THETAC 


Figure IV-2 (Cont'd) 























Figure IV-2 (Cont'd) 
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h = -ALT 

P = PATM 
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Figure IV-3. Flowchart for output subroutine 






V. TEST CASES 


In this section two test cases of the simplified computer program are 
presented. The complete input data and samples of the printout are tabulated. 
Plots of chamber pressure and thrust versus time are presented and compared 
with the results of others. Similar comparisons are made for weight calcula- 
tions and related calculated values. 

The first test case is for a 156 inch diameter SRM and the second for a 
31 inch diameter SRM. For the latter, actual static test results are avail- 
able for comparison. For the former, the comparison must be based on other 
theoretical results because the SRM has not been constructed. 

Test Case 1 - SRM Number 1 


SRM 1 is based on a design presented in Reference The basic design 
is depicted on Figure V-l and input data is tabulated in Table V-l. This 
design is of the general type under consideration for use by NASA as the 
Space Shuttle Booster. This particular SRM was selected because a relatively 
complete set of the required input data is presented in the reference as well 
as prediction results which are needed for comparison with the results of 
the simplified computer program. Also, the design is one which will test 
the use of the program for a combination circular perforated and star grain. 

The SRM has three center segments consisting of circular perforated 
grains , an aft segment with a circular perforated grain and a forward 
segment with a truncated (slotted tube) star grain. Representation of the 
grain design is a straightforward process in this case. The aft dome must, 
however, be represented by a conical closure which is specified by selection 
of the angle 0 cn . Also, since a star grain is located in the forward dome an 
effective grain length must be specified and the star grain treated as being 
flat at both ends. Since this can introduce substantial error an adjustment 
is made by introduction of tabular input values. In this case the tabular 
inputs were based on a simple linear approximation to the effect of the 
curvature using as a guide the values of star grain perimeter obtained from 
a previous run in which the effect was neglected. The effect of the cutback 
in the grain inside diameter to accommodate the submerged nozzle was neglected 
in this case but the taper in the remaining portion of the grain, which 
substantially affects tailoff characteristics, was taken into account by 
specification of and L^. 

The value of nozzle throat diameter used was 47*00 inches which differs 
from the 45.04 inches given in Reference 4 . This substitution was made based 
on information received from NASA which indicated that the higher value had 
been used in predicting the performance data given in Reference 4 . The program 
was tested with the erosion constants equal to zero because it is apparent that 
the predictions of Reference 4 were based on little or no propellant erosion. 
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A simple printout of the transients values for the initial portion of the 
trace computed by the program is given in Table V-2, and on Figure V-2 the 
chamber pressures and thrust versus time results are plotted and compared 
with results taken from Reference 4 . The reference gives only one chamber 
pressure which is presumably the nozzle end stagnation pressure. It will be 
noted that the results for nozzle end stagnation pressure and thrust are in 
agreement within 5$ except for the period between 110 to 125 seconds where 
the variation is somewhat greater. A large portion of the larger variation 
can be attributed to the approximations made in representing the effects of 
the closures. The remaining portion of the variation here as well as the 
small difference throughout the traces are more difficult to explain but are 
probably attributable to differences in analytical technique although a 
possibility exists that there could be some differences in the input data used. 

It is notable that there is good agreement in the web time calculated 
by the several approaches. Weight and related calculated values are tabulated 
as Table V-3. Comparison with corresponding data from Reference 4 are made 
later in this section. 

Test Case 2 - SRM Number 2 


SRM 2 is based on design of the TX354-5 rocket motor as presented in 
Reference 9- The basic design is depicted on Figure V-3 and input data is 
tabulated in Table V-4. This SRM was selected because of the availability 
of a relatively complete set of input data as well as specific test results. 

The grain has, in general, a cylindrical port with two radial slots. The 
main propellant grain geometry was represented by the circular perforated 
grain equations. The effect of the forward closure was represented by 
selection of 0 c h. Radii at the slots were neglected except at the aft face 
of the aft slot. Here the geometry of the entire portion of the propellant 
aft of the aft slot was represented by tabular input values because the step 
change in port area does not permit accurate representation through the 
equation inputs. In specification of x^ a the effect of the tapered insula- 
tion thickness at the bottom of the aft slot was taken into account so that xT a 
was taken somewhat larger than it would have been based on the inside taper 
alone . 

As data on the erosive burning constants of the propellant was not avail- 
able, arbitrary values were devised to give a moderate erosion rate. These 
values are listed under ALPHA and BETA in Table V-4. The tests were also run 
with ALPHA and BETA both assigned zero values; i.e., no erosive burning. 

A sample printout of the transients values calculated for the final portion of 
the trace is given in Table V-5. These are for finite values of ALPHA and BETA. 
On Figure V-4 the head end pressure and thrust versus time results are plotted 
and compared with static test results reported in Reference 9» The chamber 
pressure results of the typical test given in the reference are in this case 
presumably head end pressure since this is the pressure routinely recorded 
in standard test firings. 
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Comparison of the program output with the test data of Reference 9 shows 
that both sets of assumed values for erosive burning constants give only 
fair representation of test data for the initial portion of the trace. (It 
should be noted that the initial peaks in the test results may be at least 
partially attributable to the ignitor, effects of which are not considered in 
the program). The agreement for 'the remaining portion of the trace is 
again within 5# except for a small portion of the trace shortly before tail- 
off. The larger variation near tailoff may be traced directly to the approxi- 
mation of the effect of the forward closure made by the specification of 0 c h' 
Here an adjustment could be easily introduced to correct precisely or approxi- 
mately for the effect through use of the tabular input values . This was not 
done in the test case so that the result of the approximation could be 
evaluated by the reader. 

Again there is good agreement in the web time calculated by the program 
and that given by the reference. The agreement is poorer in case of action 
time probably due principally to neglect of the effect of chamber wall heating 
on chamber depressurization in the simplified program. 

Comparisons of Weights and Related Values 

Weights and related values for both test cases are tabulated in Table V-3. 
The values for SRM 2 are based on the finite erosion rates. 

Generally excellent agreement is found between the program and the 
reference results for propellant weights and impulse values for both SRMs. 

For SRM 1, a substantial difference in the value of P meCl p exists which results 
in rather large differences in the case thicknesses, case weights and nozzle 
weights. It is believed that this may be principally the results of differences 
in the ballistic variability factors used in the two predictions. The rather 
large difference in insulation weights for SRM 1 is attributed to neglect of 
requirements for final thermal protection in the present program. (See step 6h , 
Section III and W & , Section II). 

For SRM 2 substantial differences exist in the case and nozzle weights 
resulting in a similar difference in total motor weight. The difference in 
case weight is attributable to two factors: l) the use of approximately 10 

percent additional basic chamber thickness in the actual motor to allow for 
manufacturing variations, and 2) the use of a thick wall case segment (0.25 
inches in thickness) to accommodate a thrust ball bracket and related hardware. 
The difference in nozzle weight is principally the result of the assumption for 
the present program of a nozzle design radically different from that used in 
the TX35**-5. (See and K 2 , Section II). 

Additional testing of the program was accomplished in the process of 
compiling the design performance data given in the Appendix. 




Figure V-l. Sketch of SRM Number 1 (Reference It) 
All dimensions are in inches. 


w 



TABLE V-l. INPUT DATA FOR SRM NO. 1 


P&QP_ELLANJL.CHAEL1CJ.EEJ_S: 

RHO= 0.001990 


N = 0.250 

-^AL RHA= 0. 0_ 

BETA* 0.0 

— MU.= L.27Q0F -Q7 

CST AR= 5.16606 03 


BASIC -MQ.IQP._H IM.Q1 

L= 1243. 

— LAU= 50. A7 

VCI= 3.22000 06 


DE= 1.4B00E 02 
— DlX J_r — ft- . .7 0 . OlO.E — CU 

THETA* 0.0 


eASIC STAR GEOMETRY 

-NS= U 

LGS I = 148.00 

NP= A . 

RC= 77.470 
2.50 


NN= 0. 


TRUNCATED STAR GEOMETRY 

PP= 27.000 
T AUS = 9.67 0 


1.5000E 02 


INERT WEIGHT INPUTS 


1 0 = 0.0 




DELT AY = 0.250 
X QUT= 10 QP. O P 
DPOUT= 500.00 
7 FT A F = rt.Cpfs 


TB= 130.0 



GRAIN CONFIGURATION 

INPUT* 3 

IN* 3 

R= 2 _ 

NT= 0. 

R = 



C.P. grain geometry 

4.9 40 

5 A, 000 


S= 7. 

T NET AG= 0.0 
LGCI= 1048.50 

Jn ON I= QjlQ 

THETCN* 0.75035 
T HET C.H5._1., 5 7015. 


I DTEMP= 2.0000E Cl 

MAP* 3. 


SIGMAS* 2.0CO0E-C2 




N2= 0.0 

SYCNGM* 1 .9000 E 05 
OCC= 1.5550E 02 




DELC= 2.6200E-01 




NS EG* 5.0000E 00 
CN= 6. 2000 E 01 


SYNNOM* 1 . 90006 05 


3tI*I«M>I«l 


PSIA= 2.0000E 00 
K 1= 2 *08006-01 


K 2 - 9.2500E-02 

PS II NS = 2 .0000 E 0 


0E-02 
30Q0E-02 
T AUL= 6.50006-02 
WA= 0.0 


TABULAR VALUES 

' Y ABPK 

9,000 

9.67 9,000 

frO.O 0 

160.0 0 

All other A r s = 0 
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TABLE V-2. SAMPLE PRINTOUT OF TRANSIENT VALUES FOR SRM NO. 1 i 


INITIAL 

REYNOLDS i 

NUMBER = 

1.1206E 09 







T IME- 

0.0 

Y * 

0.0 







RN02 = 

3.9230E-01 

RHEAD= 

4.Q029E-01 

PQNOZ= 

9. 2 1 64 E 

02 

PH EAD= 

9.9906E 

02 

PTAR = 

I. 3 20 IE 00 

MNOZ = 

5.21996-01 

SUMAB = 

3. 9086 E 

05 

SG= 

4.8836E 

02 

PAT M = 

1.4696E 01 

CFVAC = 

1.7775E 00 

F VAC = 

2. 7270E 

06 

F=" 

2.4792E 

06 


TABULAR VALUES FOR VT* 9,670 READ IN 


T IME = 

0.6 

Y = 

0.250 







R NOZ = 

3.9198E-01 

RHEAD* 

3 ♦ 9945E-01 

PONOZ= 

9.1861E 

02 

PHEAD= 

9.9066E 

02 

PT AR - 

1.3434E 00 

MNO Z = 

5.0941E-01 

SUMA8 = 

3. 91 75 F 

05 

SG= 

4.8925F 

02 

P AT M= 

1.4696E <5T“ 

CF VAC - 

1.7774E 00 

F VAC = 

2. 71 98 E 

06 

F= 

2. 472 IE 

0 6“ 


T I ME = 

1.3 

Y = 

0.500 


• 




1 

K NOZ = 

3.9198E-01 

RFEAD = 

3. 989 8E-01 

PONOZ= 

9.18616 

'02 

PH 6AD= 

9.8601E 

02 i 

PT AR = 

1.3669E 00 

MNOZ = 

4.9745E-01 

SUMAB* 

3. 9263 E 

05 

SG= 

4.90 15E 

0 2 1 

PAT M- 

1.4694E 01 

C F VA C = 

1.7773E 00 

F VAC = ' 

2.72166 

06 

F- 

2.4739E 

06 i 


4 


1 T IKE = 

1.9 

Y = 

0.750 






1 

k NOZ = 

3.91986-01 

R HEA D = 

3. 9 855 E- 01 

PONOZ®*" 

■"9. 186 IE 

02 

PHEAO= 

9.8176E 

02 j 

PT AR = 

1.3906E 00 

MNO Z = 

4.860 6E-01 

SUMAB = 

3. 9351 E 

05 

SG= 

4.9107E 

02 

PAT M= 

1 .4691 E~0 1 

C F V A C = '* 

1.7 7 726 0’ 6 

F VAC = 

2. 7235 £ 

06 

F* 

2.4758E 

06 


T IME* 

2.5 

Y = 

1.000 





"TTFTOZ = 

3T5TTB E-"6Y~ 


"T796W-'0T~ 

“pimois- 

' 972053' T ~UT~ 

“phTao= 

“9rmirirr“ 

PT AR = 

1.41466 00 

MNOZ = 

4 « 7520E-01 

SUMAB= 

3.9437E 05 

SG= 

4.9199E 02 

PATM=“ 

1.4686? Cl 

CFVAC = ' 

“ 1.77 7 IE 00 

F VAC = 

2. 7310 E 06 

F= " 

2 .48 34E '0 6 


T I M E = 


Y = 

1.250 







HKT7o2>“ 

3.9218E-01 

R HEAO = 

3 • 9799E- 01 ” 

pOnoZ- 

97T55TT 

02 

PHfA^ 

“97763 IE 

ITT 

; PT AR = 

1.4387E 00 

MNO Z = 

4 • 6 52 fc£~OI 

SUMAB- 

3. 9523 E 

05 

S G= 

4.9293E 

02 

| PAT M = 

1.4 6’gOTUT 

CF VAC= 

* 

-J 

C3^ 

m 

O 

O 

F VAC - 

2. 732 BE 

06 

F= 

2.4853E 

06 







TABLE V-3. PRINTOUT OF WEIGHTS AND RELATED VALUES 


SRM Number 1 


WP1* 1 .2269E 06 

X£2* L.22 6L4£— 0..6 

WP= 1.2276E 06 

PHMA&5 9* 9SC6E_02 

ISP* 2.6215E 02 MAX EXPECTED PRESSURE* 1.1391E 03 

ISPVAC* 2.726CE 02 CYLINDRICAL CASE T HICKNESS* 6.5258E-01 

ITOT* 3.2 1 83E 08 •.{ CASE WT* 1.2395E 05 ~ 

J.TVAC* 3.34 65 E Q £ NCZZLE WT* 1.2586E 04 

I NSULAT ION W T * “~9.34G2E 03 

L INER WT* 1. 3017 E 03 

TOTAL MOTOR WT* "1.3748E 06 
ZETAM* 8.9295E -01 

RATIO OF ITOT TO WM= 2.3409E 02 ^ 

S> 

i 


SHM Number 2 


&P-.L- 3_._2.L13 e , 03 

WP2= 8.1896E 03 

WP = 8. 2035E 0 3 

PHMAX* 7.284 26 02 , 

ISP* 7.35S0E 02 f MAX EXPECTED PRESSURE* 8.1058E 02 

TSPVAC* 2 .6223E 02 ~ CYLINDRICAL CASE ThICkNESS* nWT9E-01 

I T OT = L.93 2 6£._ 0.6 CASE WT* S.73Q1E 02 

ITVAC* 2.1512E 06 NCZZLE WT= 1.1390E 02 

_ _ I NSULATION WT* I.0I05E 02 

LINER WT= 4.2155E 01 

TOTAL MOTOR WT* 9.0336E 03 

l fcT AM* “ 9'VO'eilE-Ol 

RATIO OF ITOT TO WM= 2.I394E 02 






TABLE V-U. INPUT DATA FOR SRM NO. 2 


£RCP.EL1_ANJ_.CHARACJ UUSX1C.S 



RHC= 0 .002010 

L- 0.0 610 




N= 0.27A 
jU_PHA.=.. 7.7 

[ 

TWmmuMSm^ 


BETA= 220.0 

MU= 1 .26D0F-07 


01= 7.500 

DELD 1= 0.9A0 


CST AR= 5.2110E 02 


S= 3. 

T HET AG= 0.0 


BASIC MOTOR 

n tmpw Sinw s 

L GC 1 = 173.00 ; 

LGN 1= 0.0 ! 


L = 200. 

T AU= 11.61 


T HET CN= 1.57075 j 

TF£TCH= .0 ..7 50 33 _ . _ 1 


VCI = 1.3000E OA 

..... VCF= 1.3000F 05 


1 '» 

INERT WEIGHT INPUTS 


DE= 2.2950E 01 
DT 1= 8 . 2900 E 00 


! PIPK= 9.0000E-0A 
DT£MF= 2. 20 00 E 01 


T HET A= 1.7A50E-01 
AJJLAJYr. 2.2570F-01 


SIGMAP= 3.0000E-02 
S IGMAS = 0.0 


LT AP = 1.3000E 01 

... XT= 6.0000E-01 


Nl= 3 .0000 E 00 
N2= 0.0 


2 0= 0.0 


SYCNCM= 1.5000E 05 
DCC= 3.0900E 01 


BASIC p pp FORM AN C F 

..CONSTANTS 

PSIC= 1.2000E 00 
; OELC= 2.8300E-01 


DELTAY= 0.050 
X OUT = 9999.99 


j LCC= 1.7200E 02 
1 NS EG= l.OOOOE 00 


D P OU T = 9999.99 
Z ET AF = OLt^lft 


HCN= 1. 5000 E 01 
5YNNQM = 3.5000E OA 


TB= AO .0 

_HJB= ... 0. 


PS IS= l.AOOOE 00 
PS I A= 2.0000E 00 


GAM= 1.160 
RACER= 0 .00 AO 


K 1= 2.0800E-01 

K 2= 9.2500E-02 

PS 1 1 NS = 2.0000E 00 

DEL INS= 5 • 2000E-02 


GRAIN CONFIGURATION 
LNELLI= 3 

K EH= 3.0000E-03 
K EN= 1.2000E-02 


GF A I N= 1 
S T AR= 0 


DL 1NER= 3.3000E-02 
T AUL= 6.5000E-02 


NT= 0. 

ORDER= 2 


i WA= 0.0 | 



TABULAR VALUES 


Y ABNK 

Y 

ABNK Y ABNK 


0 1586 

1.0 1562 

2.0 15^9 

3.0 1U57 

U.O 

5.0 

6.0 
7.0 

1290 8.0 7Vf 

122*1 9-0 397 

1089 9.72 0 

9 1 * 1 * 30.0 0 


All other A* s = 0 
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TABLE 7-5. SAMPLE PRINTOUT OP TRANSIENT VALUES FOR SRM NO. 2 


! T IMF* 

36.9 

Y = 11.500 







1 RNOZ * 

3.0853E-01 

RHEAD* 3.0963E-01 

PONOZ = 

6. 1 943 E 

02 

PHEAO* 

6.20276 

02 

PTAR* 

1.2804E 01 

M.VDZ* 4.651 1 F— 02 

SUMAB = 

1 . 10236 

04 

SG= 

0.0 


PAT M= 

1.46966 01 

CFVAC = 1.71846 00 

F VAC = 

5.81796 

04 

F= 

5.23726 

~04 


T IME = 

_37 *0 _ 


Y * 11.550 







RNOZ* 

3.0743E-01 

RHEAO* ‘ 3.0754E-01 

PONOZ* 

6.1150E 

02 

PHEAD* 

6.12276 

02 

PT A R* 1.26006 

01 

MNP Z = 4. 65256-02 

SUMAB* 

1.09216 

04 

SG- 

0 .0 


PATM* 

1 .4696E 

01 

CFVAC* 1.71846 00 

F VAC* 

5. 7446 E 

04 

F= 

5 • 1639 E 

04 


#»»♦*»***»***»»*****»*« 

****TAIL OFF BEGINS**** 

*»£»*<:****»****** 

| 


. TIME* 

RNCZ = 

37.2 

2.6402E- 

pi 

Y = 11.600 

RHEAO* 2.6402E-01 

PONOZ* 

3,50846 02 

PHEAD* 

r 

3.50846 02 J 

PTAR* 

1.2796E 

01 

MNDZ = 4.6540E-02 

SUMAB* 

7. 3614E 03 

SG* 

0.0 


PATM=. 1 ..463J:E .0 1 _C E V A £ * .. _L>. 7.1 B3E.....Q.Q EJ/AC= 3*?968E....Q4 E= 2..716Q E...Q4 





****BEGIN HALF 

SECOND TRACE**** 



T 1M6* 

37.4 


Y = 11.650 






RNOZ* 

PTAR* 

0.0 

1.2796E 

01 

RHEAD* 0.0 
MNOZ* 4.6540E-02 

PONOZ* 

SUMAB* 

1.7542E 02 
0. 0 

PHEAD* 

SG* 

1 .75426 02 
_Q .0 • 


PATH* 

1 .46966 

01 

CFVAC* 1.7183E 00 

F VAC* 

1.64896 04 

F* 

1.068 IE 04 


T IMF = 

37.9 


Y * 11.6 50 






! RNOZ* 

PTAR* 

0.0 

..1.27_S.6£_ 

01 

RHEAD* 0.0 
MNH 7 = 4 . 6 540F-02 

PONOZ* 
8UMA R = 

5. 7567E-01 
0.0 

PHEAD* 

5.7567E-01 
0.0 


PATM* 

1 .46966 

01 

CFVAC* 1.7184E 00 

F VAC* 

5.41156 01 

~ F=~* 

0.0 







TABLE V-6. COMPARISON OF COMPUTED WEIGHTS AND RELATED VALUES 

WITH RESULTS OF OTHERS 


SRM 1 SRM 2 


Parameter 

Present Program 

Reference 4 

Present Program 

Reference 9 

Wp *lbm 

1,227,600 

1 , 214,327 

8204 

8220 

P h max’ psia 

999 

— 

728 

721 

^sp 5 

lbf-sec/lbm 

262.2 

— 

235.6 

237.6 

I , 

sp vac* 

lbf-sec/lbm 

272.6 

270.9 

262.2 

266.1 

I T ,lbf-sec 

321,830,000 

— 

1,932,600 

1 , 953,400 

I T vac 5 
lbf-sec 

334 , 650,000 

331,300,000 

2,151,200 

2,187,000 

P ,psia 

meop 

1139 

1 ,000 

811 

799 

t ,in. 
c 

O.656 

0.570 

0.102 

0.100 

W ,lbm 
c’ 

123,950 

102,724 

573 

971 

W n , lbm 

12,580 

11,862 

113 

342 

W. ,lbm 

9 , 3^0 

11,906 

101 

143 

W^ ,lbm 

1,302 

1,278 

42 

27 

W ,lbm* 
m 5 

1 , 374 ,800 

1 , 342,097 

9033 

9703 


* No allowance is made for W . 

a 
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VI . CONCLUDING REMARKS 


The SRM design and analysis program presented in this report is based 
in part on approximate equations and methods but gives a reasonable descrip- 
tion of SRM performance. Computational times are reasonably short. The 
designer should be able to acquire quickly the skills necessary to prepare 
input data. 

Additional comparisons should be made between program and actual test 
results to establish clearly the accuracy of the program. 

Obvious improvements possible in the program fall into two categories: 
l) several modifications may be made to reduce the number of input variables 
without compromising accuracy, and 2) the method of handling the end conditions 
can be made more accurate. Most, but not necessarily all, changes in the 
latter category are toward more sophisticated programs which it was the 
objective of this analysis to avoid. 
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APPENDIX 


DESIGN DATA COMPILATION 

Here the effects of changes in various design parameters are examined 
and the results compiled in the form of traces of nozzle end stagnation 
pressure versus time* The baseline SRM used in the analysis is the same 
as SRM 1 of Section V. Except as noted below, the effect of each design 
parameter is determined separately; that is, all other input parameters 
given in Table V-l are kept the same while computations are made for 
different values of the parameter under consideration. Also, the parameters 
selected for variation have in general only minor influence on the propellant 
volume so that the designs analyzed all load approximately the same amounts 
of propellant . The stagnation pressure at the nozzle end of the grain 
was selected for plotting versus the burning time because this characteristic 
is most descriptive of SRM performance. For example, the vacuum thrust may 
be obtained by multiplying the stagnation pressure by the product of the 
vacuum thrust coefficient and throat cross-sectional area. This product 
is a constant at a given time for all the variations considered except in 
the case of changes in the radial erosion rate of the nozzle throat. Study 
of the resulting plots should provide a guide to the designer in predicting 
at least the qualitative effects of changes in the various design character- 
istics on the performance of SRMs of the general type considered. 

Radial Erosion Rate 

Figure A-l shows the effect of the radial erosion rate (E n = RADAR) on 
the stagnation pressure at the nozzle end of the grain at various times. 

As can be seen, the effect is substantial during the latter portion of 
the trace, the higher erosion rates producing significant lover stagnation 
pressure as the burning progresses. The effect on thrust, however, will be 
much less throughout the trace because of the larger throat areas which 
accompany higher erosion rates. Of somewhat more significance is the influence 
of radial erosion rate on the burning time of the SRM. As seen from the plot, 
the burning time variation is approximately four percent over the range of 
erosion rates considered. 

Initial Grain Temperature 

This parameter was investigated by adjusting the burning rate coefficient 
(a=A) according to the relationship 


a - a exp [(77 ) (l-n)ATi] 

O p K. 
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where a is the burning rate of the baseline design, which is based on an 
initial temperature of 70°F, and AT is the variation from this temperature. 
The results are shown on Figure A-2 for various initial grain temperatures. 
The generally higher pressures and reduced burning times at the higher 
initial grain temperatures are evident from the plot. 

Number of Burning Slot Faces 

The results are shown on Figure A-3. An increase in the number of 
burning slot faces (S=S) in the circular perforated grain section of the 
baseline design produces initially higher chamber pressures because of the 
added burning surface area, but the trace is more regressive because of 
the additional decrease in propellant length with distance burned. It is 
of interest that the net effect is an essentially constant burning time over 
the range of this parameter considered. 

Number of Star Points 


Influences of the number of truncated star points (np = NP) on the 
trace of nozzle end stagnation pressure versus time for a fixed length of 
the star grain in the baseline design is shown on Figure A-k . The results 
are qualitatively similar to those of the variations in the number of burning 
slot faces (S=S) in the circular perforated grain section. 

Length of Star Grain 

Influences of changes in the initial length of the star grain (Lcsi = LGSl) 
were evaluated for a fixed overall grain length (Lq s i + Lq c ^ - 1196.50 inches) 
The results plotted in Figure A-5 confirmed the expected higher initial 
pressures and greater regressivities accompanying the greater lengths of 
star grain. 

Taper 


Effects of taper were investigated by independently varying the aft end 
taper and the taper of the main portion of the grain. The former variation 
involved changes in the difference in web thickness (Xrj a = XT) over the aft 
end tapered length and the latter involved changes in the differences in web 
thickness (z Q = Z0) over the remaining portion of the grain. The results are 
shown in Figures A -6 and A-7 for the final portion of the traces which is the 
only part where a substantial effect is produced. The pressure changes for 
the initial portion of the traces were less than 1.7 percent over the range 
of variations considered, these changes being produced by changes in the 
parameter AD^ = DELDI corresponding to and equal to the changes in the 
parameter z Q = Z0. It should be noted that the effects of z Q variations axe 
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couple d with the effects of an Xp a value of 5-0 inches which is the value 
for the baseline design. Influences of the two types of changes are similar 
except that z Q increases produce much greater regressivities and burning 
times. The longer burning times result from the burning rate reductions 
which accompany the reduced pressures. It is noteworthy that if each 
center grain segment were tapered identically with a z 0 over each segment 
equal to the z Q corresponding to a continuous taper over the entire grain 
length (the present case) the results would be essentially the same, and 
therefore the alternate design could be treated by the simplified computer 
program. 

Thrust Termination Passageways 

Figure A- 8 illustrates the effects produced by use of two thrust 
terminations in conjunction with SRM 1. The passageways were assigned the 
same diameter ( I>jip - DTP) as the nozzle throat. The angle between the axis 
of the passage and the motor axis ( 6qip - THETTP) was taken as 1*5° • The 
effective web thickness (xi I , ef . f . = TAUEFF) was assigned various values. The 
discontinuities in the resultant traces occurring at y - T>p e ff are the 
result of the method used to approximate the effects of the passage geometry 
in the simplified analysis and will not occur in actual firings. To obtain 
a better approximation, the traces should be smoothed out for some distance 
ahead of the discontinuity depending on the precise geometrical situation, 
i.e., the location of the axis of the passage with respect to the head end 
of the motor. 

Erosive Burning 

Erosive burning was investigated by arbitrary variation of the constants 
(a e -b = ALPHA and £ = BETA) in the Robillard-Lenoir burning rate equation. 

The values a = 7*7 and 3 = 88.7 are the approximate values given by Peretz 
(Reference 7) Ton a composite unmetallized 78# ammonium perchlorate - 22# 
polyester propellant. Computations were made for these values and for 
changes of 50# in each variable. The results are shown in Figure A-9 along 
with the baseline results . 
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Figure A-k. Nozzle end stagnation pressure versus time for various 
numbers of star points (SRM l) 
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Figure A -6. Nozzle end stagnation pressure versus time for 

various differences in web thickness over the 

aft end tapered grain length ( SRM l). J ' a 
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Figure A-8, Nozzle end stagnation pressure versus time with 
and without thrust termination ports (SRM l). 
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